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ABSrriACT 

This document presents the software design for the Conical Scanner 
Evaluation System. The pu"pose of this system is to support the 
performance analysis of the Landsat-D Conical Scanners, which are 
infrared horizon detection attitude sensors designed for improved 
accuracy. The system consists of six functionally independent 
subsystems and five interface data bases. This document describes the 
system structure and interfaces of each of the subsystems. The 
content, format, and file structure of each of the data bases is 
specified. For each subsystem, the functional logic, the control 
parameters, the baseline structure, and each of the subroutines are 
described. The subroutine descriptions include a procedure definition 
and the input and output parameters. 
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SECTION 1 - SYSTEM OVERIVEW 


1 . 1 INTRODUCTION 

The Landsafc-D spacecraft, which is scheduled for launch in July 1982, 
will carry two Conical Scanners. These scanners will provide Earth 
direction determination data throughout the mission and attitude 
determination and control data during the initial phase of the flight. 
The attitude data obtained by the conical scanners will be used to 
control the spacecraft to a near nominal zero pitch-roll-yaw attitude. 
The conical scanners will provide data for the backup attitude 
determination and control system in case of failure of the primary 
system. The primary system will use two NASA star trackers and 
gyroscope data. The Onboard Computer (OBC) will use data from these 
sensors to compute attitudes which are expected to be accurate to a 
hundredth of degree. 

The Conical Scanner is manufactured by Ithaco, Inc. The design of the 
Conical Scanner is intended to provide improved accuracy over previous 
Ithaco horizon scanner designs. The Landsat-D mission provides an 
excellent opportunity to assess its inflight performance because of 
the availability of sensor measurements from both conical scanners as 
well as the more accurate attitude data from the OBC. An important 
part of the plan for evaluating the conical scanner involves studying 
the actual sensor measurements with respect to measurements predicted 
from the OBC :.ttitudes using comprehensive models of the conical 
scanner and the Earth. The residual errors between the actual and 
predicted measurements will be analyzed to identify sensor biases. 
Adjustments in model parameters to account for these biases will lead 
to the inflight calibration of the conical scanners. The residual 
errors remaining after inflight calibration will provide data on the 
level of accuracy obtainable from these sensors. 
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The purpose of this document is to provide a detailed design of the 
software system that will support the performance analysis of the 
Landsat-D conical scanner. Section 1 gives an overview of the system 
and its major subsystems and data bases. Section 2 defines the 
structure and contents of each data base in the system. Sections 3 
through 8 provide detailed specifications of the major subsystems, 
including baseline diagrams, subroutine descriptions, and Input/output 
requirements. 


1.2 


SYSTEM PURPOSE 


The Conical Scanner Evaluation System (CSEJ) Is a software package 
that will serve as an essential tool In the performance evaluation of 
the Ithaoo Conical Scanners onboard Landsat-D, A detailed description 
of the overall evaluation plan Is available In Reference 1. The main 
role of the CSES In the evaluation process Is summarized In the 
following steps; 

(1) Extract conical scanner data, OBC attitudes, and ephemeris from 
the spacecraft telemetry. Store the extracted data in a data 
library for later use, 

(2) Extract Earth infrared horizon radiance profiles from the Horizon 
Radiance Data Base. Integrate the radiances over the spectral 
bandpass of the conical scanner to obtain radiance profiles for 
the Landsat-D conical scanner performance analysis, 

(3) Model the responses of the conical scanner sensor optics and 
electronics to seasonal and systematic variations in the Earth 
horizon radiance, 

(11-),' Predict conical scanner measurements using reference attitudes 
(normally fcnese attitudes are from the OBC), Store the predicted 
measurements along with the corresponding observed measurements 
in a data library for later analysis. 

(5) Generate plots of the observed and predicted measurements for 
qualitative ccmparison and review. 

(6) Generate statistics on the observed and predicted measurements 
and the computed residuals. 
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1.3 


BASIC STRUCTURE 


The CSES is made up of the following subsystems and data basos, each 
of which supports a portion of the above activities: 

Subsystems 

1. Telemetry Processor 

2. Spectral Bandpass Integrator 

3. Sensor Optics and Elec*"'onics Simulator 

4. Scanner Measurement Predictor 

5. Data Plotting and Fitting Utility 

6. Scan Path Plotter 

Data Bases 


1 . Telemetry a Base 

2. Horizon ! nee Data Base 

3. Landsat-D Radiance Data Base 

4. Triggering Heights Data Base 

5. Measurements Data Base 

Figure 1-1 diagrams the basic connections between the subsystems and 
data bases. The following section gives a brief description of each 
subsystem. 










1.4 


MAJOR SUBSYSTEMS 


1.4.1 Telemetry Processor 

The spacecraft telemetry stream contains a large volume of science and 
engineering data of which only a small subset is needed for the 
conical scanner performance evaluation. The Telemetry Processor 
extracts this subset, which includes conical scanner sensor 
measurements, OBC attitudes, ephemeris, and timing information, from 
the spacecraft telemetry and stores it in the Telemetry Data Bane. 
The data can be retrieved from the data base for use in other phases 
of the evaluation process. The Telemetry Processor performs all 
necessary data conversion and validation before entering the data in 
the Telemetry Data Base. 

1.4.2 Spectral Bandpass Integrator 

The Spectral Bandpass Integrator extracts Earth horizon radiance 
profiles from an existing data base called the Horizon Radiance Data 
Base (HRDB), and processes them so that they are suitable for use in 
the Landsat-D conical scanner analysis. The HRDB contains radiance 
profiles for all latitude bands and seasons. The radiances in these 
profiles range over a set of small spectral intervals from about 8 to 
22 microns. The Spectral Bandpass Integrator integrates radiances in 
the selected profiles over the specific spectral bandpass of the 
conical scanner. These processed radiance profiles are stored in the 
Landsat-D Radiance Data Base. They are retrieved from this data base 
as needed for input bo the Sensor Optics and Electronics Simulator. 


1-6 


1.4.3 


Sensor Optics and Electronics Simulator 


The main function of the Sensor Optics and Electronics Simulator is to 
predict the conical scanner sensor responses to seasonal, systematic 
variations in the Earth horizon radiance. The predicted responses are 
stored in the Triggering Heights Data Base in the form of horizon 
triggering heights for all scanner horizon crossing positions and 
seasons. The data are retrieved from the data base as needed for 
input to the Scanner Measurement Predictor. In addition, analyses of 
the data can be performed to study the sensitivity cf the scanner to 
various parameters. 

1.4,4 Scanner Measurement Predictor 

The primary functions of the Scanner Measurement Predictor are to 
predict conical scanner measurements and generate a dataset in the 
Measurements Data Base containing the predicted measurements. The 
Scanner Measurement Predictor uses comprehensive models of the conical 
scanner and the Earth which include adjustable parameters that 
simulate all significant error sources for the scanner measurements. 
The following list summarizes the modeled error sources; 

(1) Sensor mounting alignment and scan cone angle variations 

(2) Horizon triggering height variations and Earth angular radius 
biases 

Constant measurement biases 

(4) Output voltage conversions to telemetry counts 

(5) Earth oblateness effects 

(6) Systematic horizon radiance effects 

(7) Temperature dependence 
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selectively adjusting these parameters on successive runs of the 
Scanner Measurement Predictor leads to the Inflight calibration of the 
Landsat-D conical scanners. 

The central purpose of the Scanner Measurement Predictor is to provide 
a useful, effective tool for the performance evaluation of the 
Landsat-D conical scanners. To meet this goal, its design Includes 
several input/output options. For example, while the attitude data 
needed in the model are normally the OBC attitudes from the. telemetry 
data base, there is an option to input attitude data in the form of 
constant pitch, roll, and yaw. Another option is to generate a 
dataset containing predicted measurements for one sat of model 
parameters along with predicted measurements from a different set of 
model parameters (instead of the usual predicted vs. observed 
dataset). The purpose of this feature is to readily provide 
theoretical estimates of the effects produced by variations in model 
parameters. These features and others are described in detail in 
Section 6. 

1.4,5 Data Plotting and Fitting Utility 

The Data Plotting and Fitting Utility (DPFU) produces Caloomp plots of 
selected data from the Measurements Data Base, The utility is 
designed so that any variable can be selected from a chosen dataset to 
be plotted against any other variable from the dataset. The DPFU 
provides the following options; 

o Overlaying or stacking a series of plots 
0 Labeling axes and specifying plot titles 
o Setting scales 

0 Performing a polynomial fit of the selected data 
o Performing a finite Fourier series fit of the data 
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The DPFU computes the mean of each variable and displays it on the 
plot. Typical plots include: 

0 Predicted or observed maasurements vs, time 
0 Residual errors (difference between observed and predicted) 
vs, time 

0 Various orbit and attitude parameters vs, time 

0 Correlation plots of selected variables 

^.^.6 Scan Path Plotter 

The Scan Path Plotter plots the path of the scanner field-of-vie w 
across the Earth's surface. The scan path plot helps to provide a 
general understanding of the soan geometry and is particularly useful 
for showing the scan path on the Earth when the scanner measurements 
are affected by the presence of cold clouds. The plots can be 
overlaid on GOES Earth photographs to show the scan geometry in 
relation to meteorological conditions, 

A Scan Path Plotter was developed and used in the Seasat and Magsat 
missions. The Landsat-D Scan Path Plotter is based on this existing 
utility but will be modified to accept spacecraft attitude and orbit 
parameters in a format convenient for the Landsat-D analysis. Also, , 
the utility will be improved to model the Earth oblateness in the scan 
path computation. 
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SECTION 2 . DATA BASE DESCRIPTIONS 


2.1 TELEMETRY DATA BASE (TDB) 

2.1.1 VariablffS and Formats 

The Telemetry Data Base contains the set of variables described in 
Table 2-1. 


2.1.2 File Structure 

Each dataset in the Telemetry Data Base has the following file 
char actei 'sties; 

(1) File organization; sequential 

(2) Blooksize; 3778 bytes 

(3) Record format; Fixed-length 

(4) Record length; 3778 bytes 

(5l Record structure; Table 2-2 gives a byte description of each 
record in the Telemetry Data Base. 


2-1 


Table 2-1 Telemetry Data Base Variables 


Variable 

Name 

Description 

Type 

Units 

IR0LLF(2) 

Fine Roll Error 

I»2 

Counts 

IR0LLC(2) 

Coarse Roll Error 

1*2 

Counts 

IPITF(2) 

Fine Pitch Error 

1*2 

Counts 

IP1TCC2) 

Coarse pitch Error 

1*2 

Counts 

ISIGNU2) 

Signal Status 

1*2 

Counts 

TIME 

GMT Reading at start of each 
frame in second since Sept. 1, 1957 

R»8 

Sec 

ITEMP(2) 

Scanner Temperature 

1*2 

Counts 

IB0L0(2) 

Bolometer Temperature 

1*2 

Counts 

ISENSR(2) 

Sensor Status 

1*2 

Counts 

EPARMS(2) 

Euler Parameters - OBC Reference 
Attitudes 

R»4 

N/A 

MODES (5) 

Mode flags indicating operational 
status of the spacecraft 

1*2 

Counts 

EPHPOSO) 

ECI Axis components of spacecraft 
position 

R»4 

km 

EPHVELC3) 

ECI components of spacecraft 

R»4 

km/ sec 


velocity 
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Table 2-2 Telemetry Data Base Record Structure 


Byte Contents 

Displacement 


0-1023 

TIME (1) - TIME (128) 

1024-1027 

ITEMP(1)-ITEMP(2) 

1028-1031 

IBOLO(1)-IBOLO(2) 

1032-1041 

HODES(1)-MODES(5) 

1042-1297 

IROLLFd , 1 )-IR0LLF( 128 , 1) 

1298-1553 

IROLLF( 1 ,2)-IROLLF(128,2) 

1554-1809 

IROLLCd ,1)-IROLLC( 128,1) 

1810-2065 

IROLLn{ 1 ,2)-IR0LLC(128,2) 

2066-2321 

IPITFC1 ,1)-IPITF(128,1) 

2322-2577 

IPITFd,2)-IPITF(128,2) 

2578-2833 

IPITCd ,1)-IPITCd28, 1) 

2834-3089 

IPITCd,2)-IPITC(128,2) 

3090-3345 

ISIGNLd , 1 )-ISIGNL(128 , 1 ) 

3346-3601 

ISIGNLd ,2)-ISIGNL(128,2) 

3602-3609 

ISENSR(1,1)-ISENSR(4,1) 

3610-3617 

ISENSRd,2)-ISENSR(4,2) 

3618-3681 

EPARMSd,1)-EPARMS(4,4) 

3682-3729 

EPHP0S(1,1)-EPHP0S(3,4) 

3730-3777 

EPHVELd ,1)-EPHVEL(3,4) 
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2.2 


HORIZON RADIANCE DATA BASE <HRDB) 


The Horizon Radiance Data Base HRDB contains Earth radiance data for a 
range of small apectr&l intervals around 15 microns for various 
latitudes, times of year, and Earth viewing angles. The HRDB is 
described in Reference 2. 


2.2.1 Variables and Formats 

The HRDB contains the variables described in Table 2-3. 

2.2.2 File Structure 


The HRDB has the following characteristics 


1 . 

2 . 

3. 

4. 

5. 

6 . 


7 . 


Fi7e organization ; 
Blocks! ze: 

Record format: 
Record length: 
Number of tracks; 
Number of records: 
record 1 


Direct Access 
7244 

Blocked fixed-length 
7244 
217 
217 

header record 
records 2 to 217 - data records 
Record structure; 

a. The structure for the header record is shown in Table 2-4 

b. The structure for the data records is shown in Table 2-5 


The data for a particular latitude and time of year spans two records 
on the dataset. The starting record number for a particular latitude 
and time of year (ELj^, Tj) is determined from the following equation: 

Record number =2*t(j-1) *NL+i] (2-1) 
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Table 2-3 


HRDB Variables 


Variable 

Name 

Description* 

Type 

Units 

NF 

Number of wavelengths (71) 

1*4 


NH 

Number of tangent height 
values (^11) 

1*4 


NZ 

Number of zenith angles values(IO) 

1*4 


NL 

Number of latitudes (9) 

1*4 


NT 

Number of time values (12) 

1*4 


W(i> 

(i=1,NF) 

Wavelength 

R*4 

micrcns 

H(i) 

(i=1,NH) 

Tangent height 

R*4 

km 

Z(l) 

(i=1,NZ) 

Zenith angle 

R*4 

degrees 

EL(i) 

(i=1,NL) 

Latitude 

R*4 

degrees 

T(i) 

(i=1,NT> 

Time of the year 

L*4 


R(i,k) 

Radiance 

R*4 

watt 

cm^.SR. 

(i=1,NF 





k=1 .NH+NZ) 

* The values given in parentheses are the default values. 
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Table 2-J] 


HRDB Header Record Structure 


I 


1 


Byte 

Displacement 

Name 

Description* 

0 

NF 

Number of wavelengths (71) 

J] 

NH 

Number of tangent height values (41) 

8 

NZ 

Number of zenith angle values (10) 

12 

Nl 

Number of latitudes (9) 

16 

NT 

Number of time values (12) 

20 

W2 

• 

First wavelength in microns (8.0) 

2H 

Second wavelength (8«2) 

304 

• 

''^NF 

NF-th wavelength (22,0) 

308 

»1 

H2 

• 

First tangent height in Km (0.0) 

312 

Second tangent height in Km (2.0) 

4V2 

• 

«NH 

NH-th tangent height in Km (80.0) 

476 

^2 

• 

First zenith angle in degrees (90,0) 

480 

Second zenith angle in degrees (100.0) 

516 

• 

• 

^NZ 

NZ-th zenith angle in degrees (180.0) 

520 

EL^ 

■ 

First latitude (80.0) 

556 

« 

• 

NL-th latitude (-80.0)+ 

560 

t 

First value of time (JAN.) 

608 

• 

'^NT 

NT~th value of time (DEC. > 

* Values given 

in parenthese 

are the default values. 

Negative value implies south latitude 
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Table 2-5 


HRDB Data Record Structure 


Note; Two consecutive records contain data for a particular latitude 
and time of year. 


Byte 


Displacement 

Name 

Description 


0 

RO,1> 

Radiance 

for 

(WpH,) 
(Wp.H ) 
(W3,H^) 

4 

R(2,1> 

Radiance 

for 

8 

• 

R(3,1) 

* 

Radiance 

for 

280 

R(NF,1) 

Radiance 

for 


284 

e 

na,2) 

« 

Radiance 

for 

(wJ.HgT 

• 

• 

7240 

• 

R<26,36) 

Radiance 

for 

^'^26 '*^36^ 

End of the first 

record 




BeRinning of the 

second record 



7244 

■ 

R(26,37) 

• 

Radiance 

for 


• 

• 

7424 

• 

R(NF,NH) 

Radiance 

for 

(Wujt.Hmh) 

7428 

• 

R(1,NH+1) 

k 

Radiance for 

• 

• 

14480 

• 

« 

R(NF,NH+NZ) 

Radiance for (Wjjp,Zj ^2 
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2.3 LANDSAT-D RADIANCE DATA BASE (LRDB) 

2,3,1 Variables and Formats 

The LRDB contains the variables shown in Table 2-6 


2.3.2 File Structure 

Th8 data baae «111 conalat of an unformatted direct aooeas file with 
following file oharaoteristios. 


1 . 

2 . 

3. 

4. 

5. 

6 . 

7. 


File organization: 
Blocksize: 

Record format: 
Record length: 
Number of tracks: 
Number of records: 


Direct Access 
204 bytes 
Fixed-length 
204 bytes 
4 

110 

There are two file 

2-7. 


Record structure: 
contents shown in Table 


header records with the 


The data record number for a particular latitude and time of the year 
(EL^, Tj) is determined from the following equation 

Record number s (J - 1)*NL+l+2 


The data record format is shgwn in Table 2-8. 
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Table 2-6 LRDB Variables 


Name 

Desoriptlou 

Type 

Unit 

NH 

Number of tangent heights values 



NZ 

Number of the zenith angles 



NT 

Number of times of the year 



NL 

Number of latitudes 



H(l) 
1=1, NH 

Tangent heights 

R»l| 

km 

Z(l) 

i=1,NZ 

Zenith angles 

nn 

degrees 

EL(i) 

i=1,NL 

Latitudes 

R»il 

degrees 

T(i) 

i=1,NT 

Times of the year 

I«H 

months 

RAD(i) 

1=1,(NH+NZ) 

Radiance Profile 

R»4 

watt 

cm^.SR 


Table 2-7 LRDB Header Records 


Header Record 1 


Byte 

Diaplaoement Variable 


0 

4 

8 

12 

16 

16+4*NL. 


HH 

NZ 

NL 

NT 

EL(i), i= 
T(l), 1=1 


Header Record 2 
Byte 

Diaplaoement Variable 

0 Hd), i= 

Z( 1) »i=1 


1, NL 
, NT 


, NH 
,NZ 
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Table 3-8 

LROB Data Record 


Byte 

Dlsplaocment 

Variable 

Deaoription 


0 

RAD(1) 

Radiance for H(l) 


• 

RAD(2) 

« 

Radiance for K(2) 


• 

9 

• 

» 

RAD(NH) 

Radiance for H(NH) 


4»»NH 

RAD(NH+1) 

Radiance for zenith angle 

Z(l) 

4“NH+4 

RAD(NH+2) 

• 

Radiance for zenith angle 

Z(2) 

* 

• 

4*NH+NZ-4 

• 

• 

RAD(Nri+NZ) 

Radiance for zenitli angle 

Z(NZ) 


1 



2,i| TRIGGERING HEIGHTS DATA BASE (THDB) 

2.4,1 Variables and Formats 

The Triggering Heights Data Base (THDB) contains just one variable, T 
of type R*4, which provides the horizon triggering heights in 
kilometers. 


2,4.2 File Structure 


The THDB has the following file characteristics; 


File organization; 
Blockslze; 

Record format; 
Record Length; 
Number of records: 
Number of tracks: 
Record Structure: 


Direct Access 
2896 bytes 
Fixed-length 
2896 byteT 
12 
6 

The structure of each record is shown in 
Table 2-9. 


Each record will provide the triggering height data for one month. 
The record number will correspond to the month number for which data 
is provided. 
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Table 2-9 THDB Record Structure 


Byte 

Displacement 

0 


Name Description 


T(i,j,k) Triggering heights, where l:^1,181 

indicates the satellite orbit position 
as degrees of true anomoly from the 
ascending node over the range 0 to 360 
with increments of 2 degrees. 

j = 1,2 indicates the horizon crossing, 
1 for Earth-in and 2 for Earth-out. 

k s 1,2 indicates the scanner number. 
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2.5 MEASUREMENTS DATA BASE 

2.5.1 Variableo and Formats 

The Measurements Data Base contains the set of variables described in 
Table 2-10. 

2.5.2 File Structure 

Each dataset in the Measurements Data Base has the following file 
characteristics: 

(1) File organization: Sequential 

(2) Block size: 7236 

(3) Record type: Blocked Fixed length 

(4) Record length: 108 bytes 

(5) Record structure: Table 2-11 gives a byte description of the 


record structure. 


Table 2-10 Measurements Data Base Variables 


Variable 

Name 

Description 

Type 

Units 

TIME 

Time a.')Sociated with each 

m 

Sec 

WIDTHC2) 

measurement in seconds since 
Sept. 1, 1957. 

Observed Earth width angle 

R»i| 

Deg 

PHASE (2) 

Observed Earth phase angle 

R»4 

Deg 

WIDTHP(2) 

Predicted Earth width angle 

R»4 

Deg 

PHASE? (2) 

Predicted Earth phase angle 

R»i| 

Deg 

B0L0(2) 

Bolometer temperature 

R»4 

°C 

TEMP(2) 

Scanner temperature 

R»H 

°C 

ISIGNL(2) 

Signal Status 

1*2 


ISENSR(2) 

Sensor Status 

1*2 


EILAT(2) 

Predicted Earth-in horizon 

R*4 

Deg 

EOLATC2) 

crossing latitude 
Predicted Earth-out horizon 

R*4 

Deg 

PITCH 

crossing latitude 
Reference attitude pitch 

R*4 

Deg 

ROLL 

measurement 

Reference attitude roll 

R*4 

Deg 

YAW 

measurement 
Reference attitude yaw 

R»4 

Deg 

ALT 

measurement 
Spacecraft altitude 

R*4 

km 

SUBLAT 

Subsatellite latitude 

R*4 

Deg 

SUBLON 

Subsatellite longitude 

R*4 

Deg 

ORBPHA 

Subsatellite orbit phase angle 

R»4 

Deg 

NOTE: There 

relative to the ascending node 
are two values associated with 

each variable 

that ii 


dimensioned size 2, one for each scanner. 
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Table 2-11 Measurements Data Base Record Structure 


Byte 

Displacement Content 


0-7 

TIME 

8-15 

WIDTH (1), WIDTH (2) 

16-23 

PHASE (1), PHASE (2) 

24-31 

WIDTHP(1),WIDTHP(2) 

32-39 

PHASEPd) ,PHASEP(2) 

40-47 

BOLOd) ,B0L0(2) 

48-55 

TEMPd) ,TEHP(2) 

56-59 

ISIGNLd),ISIGNL(2) 

60-63 

iSENSRd) ,ISENSR(2) 

64-71 

EILAT(1),EILAT(2) 

72-79 

EOLAT(I) ,EOLAT(2) 

80-83 

PITCH 

84-87 

ROLL 

88-91 

YAW 

92-95 

ALT 

96-99 

SUBLAT 

100-103 

SUBLON 

104-107 

ORBPHA 
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SECTION 3 - TELEMETRY PROCESSOR 
3.1 FUNCTIONAL DESCRIPTION 

The Telemetry Processor (TP) extracts attitude, ephemeris, timing, and 
conical scanner data from the spacecraft telemetry. Table 3-1 
provides a complete specification of these data by variable name, 
description, sample rate, telemetry ID, and location in the telemetry 
stream. Note that the spacecraft telemetry is structured such that a 
minor frame contains 128 8-bit words and a major frame contains 128 
minor frames. See Reference 3 for further information about the 
structure and content of the Realtime Spacecraft Telemetry (Mission 
and Engineering Formats), 

Two important functions of the TP are to attach time tags to the 
extracted data and to det^^t any anomalies in the timing information 
(e.g., unexpected gaps or backtracking). The Digital Processing Unit 
(DPU) time data (entry 15 in Table 3-1) are referenced to the major 
frame synchronization pulse that signals the beginning of each major 
frame. The TP uses this initial timing information together with the 
data rate and frame counter to time tag each row of the telemetry 
array. Specifically, at the 8K-bps data rate, each increment of the 
frame counter represents 0.125 seconds. At the IK-bps data rate, 
each increment represents 1 second. The TP uses the GMT status and 
update reports from the OBC (entry 25 in Table 3-1) as an additional 
timing and synchronization check. 

Before entering the extracted data in the Telemetry database, the TP 
performs all necessary data conversions and validations. The 
validation procedure checks that each data element is within an 
acceptable range. If a value is not within its given range, the value 
is replaced with the missing value indicator. 
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Table 3-1 Conical Scanner Evaluation 


VARIABLE 

NAI4E 


IQUAL 

IRATE 

IFHT 

f^FRA!-lE 

IROLLF(l) 

IROLLFC2) 

ISIGNL(l) 

ISIGHL(2) 

IROLLC(l) 

IROLLC(2) 

IPITF(l) 

IPITF (2) 

IPITC(l) 

IPITC (2) 

DPUTIM(6) 

ITEMP(l) 

ITEMP{2) 

IBOLOd) 

IBOLO{2) 

ISENSR(l) 

ISENSR{2) 


description 


Quality flag: 0 = good minor frame 

FF = bad minor frame 

Data rate (1 Kbps or 8 Kbps) 

Format ID (Mission or engineering) 

Frame counter 

ESA-1 Fine Foil Error 

ESA-2 Pine Roll Error 

ESA-1 Signal Status 

ESA-2 Signal Status 

ESA-1 Coarse Roll Error 

ESA-2 Coarse Roll Error 

ESA-1 Fine Pitch Error 

ESA-2 Fine Pitch Error 

ESA— 1 Coarse Pitch Error 

ESA— 2 Coarse Pitch Error 

DPU Time Code Data: Spacecraft ID, Day, Hour 

Min, sec. Msec 

ESA-1 Temperature 

ESA-2 Temperature 

ESA— 1 Bolometer Temperature 

ESA— 2 Bolometer Temperature 

ESA-1 Sensor Status 

ESA-2 Sensor Status 


Telemetry Extract (1 of 2) 


TELE- 
METRY 
JSER ID 

SAI'IPLS 

RATE 

MINOR FRAME 
FJORD 
NUMBER 

MINOR frame 
tiUMBER 
fSUBCOMS) 

C/DH 

128 

0 


C/DH 

128 

3 


C/DH 

128 

3 


C/DH 

128 

65 


ESAM 

128 

72 


ESAM 

128 

73 


ESAM 

128 

38 


ESAM 

128 

102 


ESAM 

128 

104 


ESAM 

128 

105 


ESAM 

128 

8 


ESAM 

128 

9 


ESAM 

128 

40 


ESAM 

128 

41 


DPU 

1 

32 

0-6 

ESAM 

1 

33 

44 

ESAI^ 

1 

33 

46 

ESAM 

1 

33 

47 

ESAM 

1 

33 

48 

ESAM 

4 

99 

6,38,70,102 

ESAM 

4 

99 

22,54,86,118 






Table 3~1 Conical Scanner Evaluation System Telemetry Extract (2 Of 2) 


VARIABLE 

NAME 

DESCRIPTION 

EPARHS (4) 

Euler Parameters that specify vehicle orientation 
relative to ECI frame 

MODES (4) 

Mode flags indicating stellar acquisition status, 
earth-pointing status, ephemeris source, sun 
presence 

EPHPOSO) 

ECI axis components of FS position vector 

EPHVEL{3) 

ECI axis components of FS velocity vector 

FSTIME 

Flight software time 

GMTUPD{3) 

GMT update report 




LE- 


t-lETRY 
SER ID 



IHNOR FRAME 
VJORD 
NUMBER 



115-126 9,41,73,105 

91, 95,118 23 

122 


91-95,108- 

114 

115-126 

91-95 

91-95, 

108-113 


17,49,81,113 

17,49,81,113 

30,62,94,126 

97 
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The TP generabea a report of the oharaoteristios of the output 
dataset, Including information about time span covered, the data rate, 
gaps in the telemetry stream, clusters of bad minor frames, mode 
changes, the total number of bad minor frames, and the total number 
of major frames in the dataset. 
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INPUTS AND OUTPUTS 


^ ? 


3.2.1 


NAMELIST Inputs 


The following program control parameters are Input through the 
NAMELIST input file, iTPIN: 


IFIRST 

LAST 

IHEX 


ITLM 

MAXBAD 

IDBG 


The number of the first record in the Telemetry file to 
be read. 

The number of the last record in the Telemetry file to 
be read. If LAST=0, the entire file will be read. 

If IHEXsl, each record of the Telemetry file that 
contains the OBC reference attitudes is written in the 
debug output file in hexadecimal format. If IHEX=2, 
every record of the Telemetry file is written iw the 
debug file in hexadecimal format. If IHEXsO, no debug 
output is generated. 

The logical unit number of the Telemetry file. 

The maximum number of bad records per major frame that 
are allowed before a diagnostic message is issued. 

The logical unit number of the Debug output file. 


3 . 2.2 


Dataset Inputs 


The input to the Telemetry processor is a Landsat-D Telemetry tape 
which contains the Landsat-D Realtime Spacecraft Telemetry in mission 
and/or engineering formats. The data format of this tape is defined 
in Reference 3. 
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3.2.3 


Dataset Outputs 


The telemetry processor writes to the Telemetry Data Base. See 
Section 2.1 for a description of the format and content of entries in 
the Telemetry Data Base. 

3.2.4 Printouts 


The telemetry processor generates the following printouts. 


S ummary Report 

The summary of the data span processing contains the following; 

0 NAMELIST Parameters 
0 Initial and Final GMT reading 
0 Initial Data Rate 

0 Gaps in data and changes in data rate (listed by major frame 
number, minor frame number, last GMT reading, current GMT reading, 
current data rate) 

0 Clusters (i.e., more than MAXBAD per major frame) of records that 
have been flagged as bad (listed by major frame number and minor 
frame numbers of bad records), 
o Total number of major frames in output dataset 
0 Total number of bad records 

0 Mode flags and indication of any changes in the mode flags 
Debug Output File (Optional) 


Contains selected telemetry records of the input and output datasets 
written in hexadecimal format and selected internal variables. The 
amount of debug output is controlled by NAMELIST parameters. 
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3.3 


SYSTEM STRUCTURE 


Figure 3-1 aontaina fche baaeline diagram of the Telemetry Prooeaaor. 
The routine TPMAIN aervea aa the driver of the ayatem. It initializea 
all neoeaaary paramotera, performa all input and output, and calla 
UNPACK, TIMTAG, and CHECK to parae the telemetry atream, perform time- 
tagging, oheok for gapa in the telemetry, and flag miaalng or out-of- 
range valuea. Figure 3-2 providea a flow diagram of the ayatem. 








ORIGINAL PACK lU 



AWOTIIER^ 
MINOR 
FRAME ^ 


CALL TIHETAO 
TO TAO EACH 
MINOR FRAME 
AND PERFORM 
TIHINC CliCCKS 


IF NEEDED, 
WRITE TIMINC 
INFO. IN j 
SUMMARY / 

REPORT / 


CALL CHECK 
TO TEST FOR 
OUT-OF-KANCE 
VALUES 


WRITE MAJOR 
FRAME TO 
OUTPUT J 

DATASET / 


another^ 

MAJOR 
FRAME , 
S.7 / 


WRITE 

SUMMARY 

REPORT 


Figure 3-2 Telemetry Processor Flow Diagram 


U - 







SUBROUTINE DESCRIPTIONS 




The purpose, procedure, inputs, and outputs of the telemetry processor 
subroutines are described as follows: 

TPMAIN 


Purpose; TPMAIN serves as the driver of the Telemetry Processor. 
Procedure: 

(1) Intializes parameters and reads NAMELIST input file, 

(2) Reads telemetry records into buffer, 

(3) If requested, writes contents of buffer in debug output 
file in hexadecimal format. 

Calls UNPACK to unpack the buffer and load data into 
appropriate arrays. 

(5) Calls TIMTAG to tag each record and perform time checking. 

(6) Calls CHECK to test for out-of-range values and flag any 
that are found. 

(7) Writes telemetry dataset in Telemetry Data Base. 

(8) Writes summary report. 

Input Parameters; None 
Output Parameters: None 
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UNPACK 


Purpose: UNPACK extracts the necessary data from the input buffer and 

stores them in the appropriate arrays. 


Procedure: 

(1) Checks the quality flag located in the first byte of the 
telemetry record buffer, BUFFER(I), If the flag is nonzei'o, 
the record is bad and missing value indicators are loaded 
into the data arrays. 

(2) Unpacks the data that are in fixed locations on every 
record, i.e,, data rate, sensor measurements. 

(3) Depending on minor frame number, unpacks selected 
subcommutator data, i.e., OBC attitudes, ephemeris, other 
OBC reports. 

Input Parameters: 

BUFFER Buffer containing one minor frame (128 bytes) of the 
telemetry stream. 

NMINOR Minor frame number. 

Output Parameters: 

IQUAL Array containing qualilty flag of each minor frame: 
IQUAL(I)=0 Minor frame I contains valid data. 

= 1 Mi.::><” frame I is missing. 

IRATE Array containing data rate indicator of each minor frame: 
IRATE(I)=1 Data rate is IK-bps 
=8 Data rate is 8K-bps 
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ICSDAT 


ITEMP 

IBOLO 

ISENSR 

EPARMS 

MODES 


128 by 10 array containing sensor measurements from 
conical scanners 1 and 2 where 
Column 1 = Scannerl Signal Status 

2 s Scanner2 Signal Status 

3 = Scannerl Roll Fine Error 

4 = Scanner2 Roll Fine Error 

5 = Scannerl Roll Coarse Error 

6 = Scanner2 Roll Coarse Error 

7 = Scannerl Pitch Fine Error 

8 = Scanner2 Pitch Fine Error 

9 s Scannerl Pitch Coarse Error 
10 = Scanner2 Pitch Coarse Error 

Array containing scanner temperatures for scanners 1 and 2 
Array containing bolometer temperatures for scanners 
1 and 2 

Array containing sensor status flags for scanners 
1 and 2 

Array containing the OBC reference attitudes which are 

expressed as Euler parameters 

Array containing mode flags from an OBC Report: 

MODES(I) s 1 Coarse sun acquisition 
= 2 Fine sun acquisition 
s 3 Stellar acquisition 
= 4 Earth pointing mode 
= 5 Inertial hold mode 
= 6 Slew mode 
= 7 Orbit 

M0DES(2) = 1 Earth pointing mode 

= 2 Stellar pointing mode 
M0DES(3) = 0 EphemerieS is from GPS 
= 1 Ephemeris is uplinked. 
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MODES(ij) = 1 Stellar acquisition in progress 
s 2 Stellar acquisition completed 

EPHPOS Array containing ECI axis components of spacecraft 

position. 

EPHVEL Array containiing ECI axis components of spacecraft 

velocity. 

NOTE: The contents of EPHPOS and EPHVEL depend on the 

value of MODESO). 

DPUTIM Array containing the DPU time which is referenced to the 
start of each major frame. 


TIMTAG 


Purpose: TIMTAG attaches time tags to each minor frame and checks for 

anomalies in the timing information. 


Procedure; 

(1) Sets the current major frame start time using the DPU time 
data which is referenced to the major frame sync pulse, 

(2) Checkd the current major frame time against the previous 
reading to detect timing gaps or backtracking. 

(3) Uses data rate information to compute time tags for each minor 
frame. 

(4) Flags changes in data rate. 

Input Parameters : 

DPUTIM Array containing DPU time daca in days, hours, minutes, 
seconds, milliseconds 

LASTIM Initial time of previous major frame. 

LSTRAT Final data rate of previous major frame. 

IRATK Array containing data rate information for each minor 
frame in current major frame. 

Output Parameters: 

CURTIM Initial time of current major frame. 

TIME Array containing time tages (in seconds since Sept. 1, 
1957) for each minor frame. 


CHECK 


Purpose: Validate data and flag missing or out-of-range values. 

Procedure: 

(1) Check quality flags. If the total number of bad records 
exceeds the specified limit per major frame, write message 
in summary report. 

(2) On each valid minor frame, test the contents of selected 
variables against validation criteria. 

(3) Update running totals of missing or invalid values for each 
variable. 

Input Parameters: 

NMAJOR Number of current major frame. 

NMISS Array containing running total of number of missing or 
out-of-range values for each variable that must be 
validated. 

IQUAL Array containing quality flag of each minor frame. 

ICSDAT 128 by 10 array containing sensor measurements from 
scanners 1 and 2. 

ITEHP Array containing scanner temperatures for scanners 1 and 

2 . 

IBOLO Array containing bolometer temperatures for scanner 1 

r id 2. 

ISENSR Array containing sensor status flags. 

EPARMS Array containing OBC reference attitudes. 

EPHPOS Array containing ECI axis components of spacecraft 
position . 

EPHVEL Array containg ECI axis components of spacecraft 

velocity. 


/; 
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Output Parameters; 


Missing value indicators are installed in the above arrays where 
needed. Updates are made to NMISS as needed. 
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SECTION 4 - SPECTRAL BANDPASS INTEGRATOR 


4.1 FUNCTIONAL DESCRIPTION 

The Horizon Radiance Data Base (HRDB) contains "average" Earth horizon 
radiance profiles. The radiances in these profiles are a function of 
four parameters: tangent height or zenith angle, latitude, time of 

year, and spectral band, i.e., 

R = F(H or Z, EL, T, W) 

where 

H = tangent height in kilometers for lines of sight that pass above 
the horizon 

Z = zenith angle in degrees for lines of night that intersect the 
Earth 

EL s Earth latitude 
T = time of year 

W = wavelength interval in microns. 

The purpose of the Spectral Bandpass Integrator is to generate a new 
data base (LRDB) containing Landsafc-D radiance profiles by integrating 
the HRDB radiances, with respect to the spectral transmission function 
of the Landsat-D conical scanner which is shown in Figure 4-1. 

The computational procedure used in the Spectral Bandpass Integrator 
involves the following steps: 

1- Using the conical scanner spectral transmission function, 
determine the weight to be applied to the radiance data for each 
spectral interval in the HRDB. The weight is the product of the 
spectral transmission and the band width i.e.. 
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Figure 4-1 Spectral Transmission Function of the Landsat-D Conical 
Scanner. 
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WEIGHT(l) = BANFUN(i)»[W(i+1)-U(i)] 

Where BANFUN(i) is the average transmission in the interval 
between W(i) and W(i>»-1), 


For each tangent height or zenith angle, latitude, and aeaaon, 

compute the weighted aum of the radianoea over the apeclral banda, 
i • e e • 

NF 

RAD(H or Z,EL,T) = z F(H or Z,EL,I,W(i))«WEIGHTa> (i^- 3 ) 

i=1 


for all H or Z, EL, and T 

where NF Is the number of frequencies or wavelengths in the HRDB 


Enter the computed radiances and corresponding parameters in the 
Landsat-D Radiance Data Base. 


4.2 


INPUTS AND OUTPUTS 


4,2.1 

Namelist Inputs 



The following liato the namollat Input paramotera Included In 

the 

NAMELIST 

&SPBAIN. 



BANFUN(i) 

(1=1,71) 

WSTART 

Array containing the average spectral 
corresponding to each wavelength Interval in 
Starting wavelength for bandpass integration 

transmission 
the HRDB 

WEND 

(for Landsat-D WSTART = 13.6 microns) 
Ending wavelength for bandpass integration 
(for Landsat-D WEND =16.4 microns) 



NHRDB 

Unit number for HRDB 



NLRDB 

Unit number for LRDB 


t 

NPRINT 

Unit number for printout 



4.2.2 

Dataset Inputs 



The input radiance data to the Spectral Bandpaas Integrator la from 
the Horizon Radlanee Data Baae (HRDB). This data baae la daaoribed In 

Section 

2.2. 




4.2.3 Dataset Outputs 

Th. output rudianoe data for tha apaolflad apactral bandpaas ta 
written to the Landaat Radiahoa Data Baae (LRDB). See Section 2.3 for 
a deaorlptlon of the format and content of the LRDB. 
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4.3 


SYSTEM STRUCTURE 


The data flow for the SPBAIN is shown in Figure 4-2. Figure 4-3 
provides a flow diagram of the system. 

4.4 SUBROUTINE DESCRIPTIONS 

The purpose, procedure, inputs and outputs of the SPBAIN subroutines 
are described as follows: 

SPMAIN 


Purpose: SPMAIN serves as the driver of SPBAIN. 

Procedure: 

(1) Read the input data from the NAMELIST and HRDB header record 

(2) Read HRDB data records for a particular latitude and Mme 

(3) Call TRAP tc provide the integration over the specified 
wavelength range 

(4) Write one LRD3 data record for a particular latitude and 
time and generate printout of the results 

(5) Repeat steps 2 through M for all latitudes and times 

(6) Write the LRDB header record 

Input Parameters: NONE 


Output Parameters: 


NONE 



















TRAP 


Purpose; TRAP is a subroutine that integrates over the specified 
range of the wavelength, 

Propfldure; Sum over the first element of a two dimensional array 
using an input array to weight each element, i.e. 


NF 

RAD(k) = 2 R(i,k)«WEIGHT(i) 
i = 1 

for all k. 

Input Parameters: 


BANFUN(i) 

(i=1,NF) 

R(i,k) 

(i=1,NF) 

(k=1, NP) 

NF 

NP 


Array containing the spectral bandpass function 
Array of radiance profiles from one HRDB record. 


Number of wavelengths or frequencies 

Total number of tangent heights and zenith angles in 

the radiance profiles. 


Output Parameters; 

RAD(k) Radiance profile integrated over the spectral bandpass 

(k=1,NP) 
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SECTION 5 - SENSOR OPTICS AMD ELECTRONICS SIMULATOR (SOES) 


5.1 FUNCTIONAL DESCRIPTION 

The Sensor Optics and Electronics Simulator (SOES) is designed to 
estimate the effects of seasonal systematic radiance variation on the 
conical scanner measurements. A key input to this program is a data 
base describing the Earth radiance as seen by the Landsat-D conical 
scanner spectral bandpass. Other NAMELIST inputs describe t)ie sensor 
optics and electronics and other modeling parameters. The key output 
is a database describing the apparent horizon triggering heights as a 
function of spacecraft position and season. This program is designed 
so that it may also be used as an anaiyt;ical bool for studying the 
effects of various detailed model parameters on the sensor 
measurements. 

A flow diagram of SOES is given in Figure 5-1. An outline of the 
computational steps is given below. 

Step 1 

The NAMELIST input is read and initialization computations are 
performed, including computation of impulse response function 
parameters (used in the convolution integral, step 5) from the 
transfer function parameters. 

Step 2 

The Landsat-D radiance profiles for a particular month are read from 
the LRDB and interpolated to 1 degree latitude increments. 
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Figure 5-1 SOES Simulation Procedure 
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step 3 


The radianoe Is integrated over the bolometer field-of-view (FOV) to 
obtain the total incoming radiance to the bolometer. The radiance is 
assumed to vary in the direction normal to the Earth horizon. The FOV 
is assumed to be shaped like a parallelogram with sharp edges and with 
form and orientation dependent on the prism rotation angle. 

Step ^ 

Steps 1 and 2 are repeated for each scan position in order to obtain 
the predicted input voltage signal as a function of time as the 
scanner sweeps across the Earth, 

Step 5 


This input voltage signal is convolved with the system impulse 
response function to obtain the predicted output voltage signal at 
representative times as the scanner sweeps across the Earth. The 
output signal will be similar to that shown in Figure 5-2. 

Step 6 

The electronincs output voltage signal as the scanner sweeps across 
the Earth is examined to obtain the maximum voltage, as the 
scanner enters the Earth and the minimum voltage, as the scanner 
leaves the Earth. Parobolic interpolation is used to obtain these 
values. These peak voltages are used to set the threshold levels for 
the next step. 
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step 7 


The threshold crossing times of the output signal, and Tg for the 
Earth-in pulse, and and Tjj for the Earth-out pulse, are determined. 

Polynomial interpolation is done on the output signal and a Newton- 

« 

Raphson iteration used to find the threshold crossing times. 

Step 8 

Based on the threshold crossing times, the scanner Earth-width and 
Earth-phase measurements are computed. These measurements are used to 
compute the apparent triggering heights at the horizon crossings. 

Step 9 

Steps 3 through 9 are repeated for each scanner. 

Step 10 

Steps 3 through 8 are repeated for each spacecraft position around the 
orbit. The results are written to the Triggering Heights Data Base 
(THDB) . 

Step 11 

Steps 2 through 10 are repeated for each month of the year for which 
radiance data is available. 
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INPUTS AND OUTPUTS 


<},2 

5.2.1 Namelist Inputs 

The NAMELIST input for SOES, &SQESIN contains the following variables. 
General Control Parameters 


IMODE Program flow mode 

= 0 analysis mode, run for a selected month and time span 
s 1 THDB generation mode, run for all twelve months for full 
oribts. lORBIT is automatically set to 1 to calculate 
orbit positions appropriate for the THDB 
IMONTH Selected month number (1 to 12) for radiance data, used when 
IMODE = 0 


Earth Modeling Parameters 
RO Nominal Earth radius 

HT Nominal horizon triggering height for ideal Earth 

Satellite Orbit Parameters 


lORBIT Orbit position selection option 

= 0 use orbital elements and TSTART, TSTEP, TSTOP 
= 1 use circular orbit with ASEMI and EYE and calculate 
orbit position appropriate for output to THDB 
TSTART Starting time (YYMMDD.HHMMSS) 

TEPOCH Orbit epoch time (YYMMDD.HHMMSS) 

ASEMI Semimajor axis (kilometers) 

E Eccentricity 

EYE Inclination 
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EHO Mean anomaly 

WO Argument of perigee 

RANODE Right ascension of the ascending node 

TSTEP Time step in calculating positions around the orbit (seconds) 

TSTOP Stop time if full orbit is not desired (seconds) 

Scanner Parameters 


AZI(l) Azimuth angles for scanners and 2 

( 1 = 1 , 2 ) 

TILT(i) Tilt angles for scanners 1 and 2 
(i=1,2) 

TWIST(i) Twist angles for scanners 1 and 2 
(i=1,2) 

GAM(i) Cone angles of FOV center from scan-axis for scanners 1 and 2 
(i=1,2) 

SPIN(i) Spin rates for scanners 1 and 2 
(i=1,2) 

Radiance Profiles Expansion Parameter 


NPTSL Number of points to be used in the Eul er-Lagr ange 

interpolation over latitude at fixed values of tangent 
height and zenith angle. 

Incoming Signal Computation Parameters 

NPHI Number of scanner rotation angles for which the incoming 

radiance is calculated (NPHI £ 200) 

PHI(i) Scanner rotation angles for which the total incoming 

(i=1,200) radiance is calculated. 


Fleld-of-Vlew Integration Parameters 


FOVSIZ Width of scanner square field-of-view prior to distortion by 
prism (degrees) 

FOVSTR Stretching factor of field-of-view due to prism distortion 
FOVOHI Rotational orientation of the square bolometer flalte with 
respect to the scanner y-z axis (degrees) 

NPTSZ Number of interpolation points for radiances to be used in 
Euler-Lagrange interpolator at fixed value of latitude and 
input value of zenith angle 

NPTSH Number of interpolation points for radiance to be used in 
Euler-Lagrange interpolator at fixed value of latitude and 
input value of tangent height 

Electronic Transfer Function Parameters 

RADVOL Radiance to volts conversion factor 

ATB Magnitude of scan angle whose negative corresponds to the 

time at the limit of the convolution integral 
TMAXIN Maximum time period for the convolution integral 
THRESH Ratio of triggering signal to maximum signal 

NPTSI Number of points to be used in the interpolation on the 

electronics input signal 
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Threshold Caloulablon Control Parameters 


El 
E2 
B1 
B2 
A1 
A2 
A3 
Ai< 

A5 
A6 

ESP Constant specifying the degree of convergence in the Newton- 

Raphson procedure for locating the threshold 
lEND Maximum number of iterations in RTNI ' 

DELTA! Time offset for the computation of the numerical derivative 
for the Newton-Raphson procedure 

NPTSO Number of points to be used in the interpolation on the 
electronics output signal. 

Input/Output Control Parameters 

IPLOT Plot input/output signal 
= 0 no plot 

= 1 plot input and output impulse 
= 2 above and intermediate result 
IFTOUT Printout unit number for result 

IFTDEB Printout unit number for debug 

ILRDB Fortran unit number for lRDB input data file 

ITHDB Fortran unit number for THDB output data file 


Coefficients of the transfer function 
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ILEVEL(i) 

(i=1,30) 


Debug level for 1-th subroutine printout 
= 0 no printout at all 
= 10 start/end messages and parameters 
= 20 above and input/ouput parameters 
= 30 above and intermediate result 
s 40 above and debug type 1 messages 
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5.3 


SYSTEM STRUCTURE 


The data flow for the SOES Is illustrated in Figure 5-3. Input to 
SOES consists of the NAMELIST parameters and Landsat-D radiance data 
(LRDB). Output from SOES includes triggering heights data versus 
subsatellite latitude which is written to the Triggering Heights Data 
Base (THDB). The SOES consists of a main program and orbit generator, 
an interpolator, an integrator, and various subprograms. The 
analytical basis for various subprogram computations is described in 
Reference 4. Figure 5-4 illustrates the program structure. 
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Figure 5- 3 SOES Data Flow 












































SUBROUTINE DESCRIPTIONS 


5.i| 


The Sensor Optics and Electronics Simulator subroutines are described 
as follows. Brief descriptions of the purpose of the existing 
subroutines are included. The existing subroutines are taken from 
Reference 5. 


SEMAIN 


Purpose: SEMAIN is the driver of the Sensor Optics and Electronics 

simulator. 


Procedure: 

(1) A card image copy of the NAMELIST input dataset is written out 
and the NAMELIST input is read. If the program is in analysis 
mode (IMODEsO) it is set up to run for a particular time of year, 
otherwise (IM0DE=1) is set up to loop though all the months of 
the year and write to the Triggering Heights Data Base. 

(2) Call subroutine RESDUE to intialize the impulse response function 
parameters. 

(3) Read the LRDB for a particular month. 

^4) Call subroutine ORBGEN to compute the spacecraft position. 

(5) Compute the scan axis attitude. 

(6) Compute the field-of-view center line-of-sight. 

(7) Call subroutine AVRAD to calculate, the radiance over the field- 
of-view. 

(8) Repeat steps 6 through 7 for the next scanner rotation position 
until the full Earth scan is done. 

(9) Call subroutine LOCATE to compute the output signal from the 
electronics, and apply the horizon locator logic to find the 
Ea^th-in and Earth-out crossing angles, 

(in) Calculate the horizon triggering height. 

(11) Repeat steps 5 through 10 for the second scanner. 

(12) Repeat steps 4 through 11 for the next spacecraft position until 
a full orbit is done or the ?■>. time is reached. 

(13) Store the results (trig ,r'" ■*' ghts) in THDB. 

(14) If IM0DE=1, repeat ste.o£ ^ ^nrough 13 for the next month until 
the full year is done. 
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Input Parameters; None 


Output Parameters?, None 


RESDUE 


Purpose: Subroutine RESDUE initialize the constants for the impulse 

response function. 

Procedure: The value of various parameters used in the transfer 

function (see formulas in Appendix A, Reference ‘D are 
passed to RESDUE via common block /POLES/ to calculate 
coefficients for the impuse response function. 

Input Parameters; 

Through common block /POLES/ - Coefficients for the 
transfer function 

El, E2, Bi, B2, A1, A2, A3, A4, A5, Au 

Output Parameters: 

Through common block /POLES/ - Coefficients for the impulse response 
function 

All, AI2, AI3, AI^l, AI5, AI6, ALFA?, ALFA8, BI1, BI2 CI1, CIH 


XSPAN 


Purpose: 

XSPAN accepts radiance data for latitudes of -80 degrees to +80 
degrees at 20 degree increments. It interpolates radiance data 
for latitudes from -90 degrees to +90 degrees at 1-degree 
increments and stores the results in 181-oolumn arrays. An 
interpolation routine AIR that uses Aitken's recursive method is 
used to perform the interpolations. 


AIR 

Purpose : 

AIR (Aitken's interpolation routine) interpolates on input "X" 
and "Y" tables by fitting a polynomial to an input number of 
points in the tables using Aitken's recursive method. The points 
chosen for the polynomial fit are those nearest the X value which 
is input for interpolation. The output is the interpolated Y 
value which corresponds to the input X value. 

XSPAN Z 


Purpose ; 

XSPANZ stores the interpolated radiance data versus zenith angle 
received from XSPAN into common blocks for use by the subroutines 
RZNPFL and RZSPFL. The data are stored in 1-degree increments 
from 90°S to 90°N for all zenith angles. 


XSPANH 


Purpose ; 

XSPANH performs the same function as XSPANZ except that it stores 
the data for the tangent heights into common blocks for use by 
the subroutines flHNPFL and RHSPFL. 

TCON40 


Purpose: 

TCONilO converts time in YYMMDD.HHHMSSMMM to time in seconds from 
September 1, 1957t at 0 hours UT, 


JD 


Purpose : 

JD converts calendar days in the format YYYYMMDD to Julian day 
(days since January 2, ^<713 BC) . 

ORBGEN 


Purpose : 

ORBGEN is a Keplerian orbit generator. It accepts time and 
orbital elements and generates the spacecraft positon and 
velocity in geocentric inertial coordinates. 
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RMAT 


Purpose : 

RMAT computes the inertial-to-orbital coordinates transformation 
matrix. It accepts spacecraft position and velocity vectors and 
outputs rotation matrices from orbital to geocentric inertial 
coordinates (GCl) and from GCI to orbital coordinates. 


UNVEC 


Purpose : 

UNVEC unitizes a vector. It accepts a vector of nonzero 
magnitude and provides a unit vector in the direction of the 
original vector and a scalar that equals its magnitude. 

SCNLAT 


Purpose : 

SCNLAT calculates the latitude corresponding to the center of the 
scanner FOV, using the scanner ID and spacecraft ephemeris and 
assuming a spherical Earth with pitch, roll, and yaw equal to 
zero . 


GHPRD 


Purpose : 

GMPRD computes a general matrix product. 
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AVRAD 


Purpose: AVRAD computes the average radiance over the scanner's FOV 

taking into account the variation in radiance over the FOV 
(Equation 4-9 through 4-10 in Reference 4), The atmosphere 
is assumed to be composed of shells of equal radiance such 
that an Isoradiance contour that is concentric with the 
Earth's surface can be defined. This reduces the field-of- 
view averaging to a single-dimension integral, perpendicular 
to the isoradianoe contour. The equations for this 
computation are discussed in Reference 4, Section 4,3.2. 


Procedure: 

(1) Call INTLIM to compute parameters U and KCI that define upper and 
lower limits for the FOV integrals. 

(2) Call SIMPSN three times to calculate three integrals with 
different limits and the different external functions FI, F2, and 
F3, to complete the field-of-view integration in three pieces. 

(3) Sum up the total radiance over the field-of-view. 

Input Parameters: 

PHIIN Scanner rotation angles 

ID Scanner identification (1 or 2) 

Through common block /FOVPAM/ 

F0VS1Z \ 

FOVSTR I 

FOVORI \ Field-of-view integration parameters 

NPTSL 1 
NPTSH 1 
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Output parameters; 


RAD 


The radiance averaged over FOV 
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INTLIH 


Purpose; Subroutine INTLIM calculates the parameters U ana KCI that 
are used by Simpson's integration procedure as upper and 
lower limits. 


Procedure; 

(1) Calculate the angle ALFA between the direction of the scan path 
and isoradiance normal vector R (Equation ^l-9t Reference *0 

(2) Calculate U and KCI (Equation 4-10 in Reference 4) 

Input Parameters; 

GAMMA Cone angle of FOV center from scan-axis 

RO Nominal Earth radius 

FOVSIZ Width of scanner square field-of-vlew prior to distortion by 
prism (degrees) 

FOVSTR Stretching factor of fleld-of-view due to prism distortion 

FOVORI Rotational orientation of the square bolometer flake with 
respect to the scanner y-z axis (degrees) 

ETA Separation between Earth center and scan axis (degrees) 

Output Parameters: 

U Distance from the FOV center to the farthest edge of the 

FOV, projected on the direction that the radiance changes. 

KCI Distance from the FOV center to the FOV corner point 

projected on the direction that the radiance changes. 
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FI. F2, F3 


Purpose: 

FI, F2, F3, are the external functions that provides the first, 
and third integral in the subroutine AVRAD in Section 5.4,13. 
Fj^ = (U+X)R(X) 

?2 = R<x) 

F 3 = (U-X) R(X) 

Input parameters: 

X “ integration variables 

U “ integration parameters from subroutine AVRAD 
Output parameters: 

F1,F2,F3 Function value evaluated at X 


second 
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SIMPSN 


Purpose ; 

SIMPSN is a utility subroutine that performs an integration over 
an input function over a specified range. 

Method ; 

The integral is evaluated using Simpson’s integration rule. 

Input Parameters; 

AA - uppper limit of integration 
BB - lower limit of integration 
NPT- number of the integration points 

FX - external function name that evaluatesthe integrand 
Output parameters; 


AHS - value of integral 


RZSPFL 


RZSPFL accepts equivalent scanner latitude and zenith angles. It 
accepts 90 radiance profile vectors covering the latitudes of 90 
degrees to 1 degree at 1-degree intervals. RZSPFL finds the 
vector closest to the scanner latitude and interpolates on the 
vector elements to find the radiance corresponding to the given 
zenith angle. It is used only for southern latitudes. 

RHSPFL 

RHSPFL is the same as RZSPFL, except that it operates on tangent 
heights rather than zenith angles, 

RZMPFL 


RZNPFL is the same as RZSPFL, but for northern latitudes. There 
are 91 columns (the Equator is included). 

RHNPFL 

RHNPFL is the same as RHSPFL, but for northern latitudes. There 
are 91 columns (the Equator is included). 


GRAPH 


GRAPH is a general line printer plot subroutine. It is used by 
SOES to make plots of the electronics input and output signals 
and the triggering heights output. 
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LOCATE 


Purpose: 

The subroutine LOCATE obtains the output signal from the sensor 
electronics, and determines the horizon crossing positions that 
are measured by the sensor. 

Procedure: 

(1) Colls the subroutine SIMPSN with the external function CONVOL to 
convolve the Input signal to the bolometer with the Impulse 
response function to obtain the output voltage from the 
electronics. The call to SIMPSN Is repeated to obtain a table of 
the output voltage as a function of time, 

(2) The subroutine SIGMAX Is called to obtain the peak signal In the 
Earth-in horizon crossing. 

(3) The voltage threshold level is determined as a fixed percentage 
of the peak signal. 

(t|) The subroutine GUESS is called to obtain initial estimates for 
the threshold crossing times. 

(5) The subroutine RTNI is called bo determine each of the two 
threshold voltage crossing times for the Earth-in crossing. 

(6) Steps 2 through 5 are repeated except to find the threshold 
voltage crossing times for the Earth-out horizon crossing. 

(7) The threshold voltage crossing times are used to compute the 
scanner rotation angles at the horizon crossings. 

Input Parameters: 

PHI Array containing the scanner rotation angles at which the 

incoming radiance is computed. 
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RADAV Array containing the radiance ,Y''eraged over FOV 
NPHI Number of points in arrays PHI and RADAV 

Output Parameters: 

PHITRG Array containing the triggering rotation angles 
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Purpose: CONVOL is a FORTRAN FUNCTION subprogram that is used by 

subroutine SIMPSN as an external function and provides the 
intregand of the convolution integral. 


Procedure: 

(1) Calls the subroutine SINPT to obtain the input electronic signal 
function at the time TPRIHE 

(2) Calls the subroutine EVAL to evaluate the impulse response 
function at the time difference between TPRIME and T which is the 
time at which the output voltage is being computed. 

(3) Calculates the product between the input electronic signal 
function for the given time and the impulse response function for 
the time shift from T. 

Input Parameters: 

TPRIME Argument of the function 
Through common block /LIMIT/ 

T Time for wh'iv.h the output voltage is being computed 


Output Parameters: None 


SINPT 


Purpose; SINPT evaluates the input voltage to the signal processing 
electronics at a given time. 


Procedure; 

(1) The scanner rotation angle, PHIT, at the input time is calculated 
using the scanner spin rate, 

(2) The subroutine AIR is called to obtain the average incoming 
radiance at the rotation angle PHIT by interpolating on the 
radiance vs. scanner rotation angle table, 

(3) The incoming radiance is converted to a voltage by using the 
NAMELIST parameter RADVOL 

Input Parameters; 

TIME Time at which to evaluate the input voltage to the signal 

processing electronics 

RATE Spin rate of the scanner 

RADVOL Radiance to volts conversion factor 

Output Parameters; 


SIGNAL 


Input signal to the electronics in volts. 


EVAL 


Purpose: Subroutine EVAL evaluates the impulse response function for 

a given time. 


Procedure: The coefficients for the impulse response function were 

calculated by subroutine RESDUE and saved in common block 
/POLES/. These coefficients are used to calculate the 
impulse response function (Appendix A, Reference <0 . 

Input Parameters: 

TPRIHE Argument of the impulse response function 

Through common block /POLES/ - Coefficients for impulse response 
function. 

All, AI2, AI3» AI1|, AI5, AI6, ALFA?, ALFA8, BI1, BI2, CI1, CI2 • 

Output Parameters: 

TRF Impulse response function calculated for the given time 
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SIGMAX 


Purpose: The subroutine SIGMAX calculates the maximum or minimum of 

the electronic output vo’,i'<age signal and obtains the voltage 
threshold level. The ar» ay containing the output voltages are 
searched until three points around the maximum (or minimum) is 
found. Then parobolic interpolation is used to find the maximum 
(or minimum) . 

Procedure; 

(1) The search procedure starts with the array index I set equal to 
the input variable INDEX1, 

(2) The three consecutive points ANS(I), AHS(I+1) AND ANS(I+2) are 
evaluated for the maximum or minimum condition. 

(3) The array index I is incremented by one and step 2 is repeated 
until the three maximum or minimum points are found. 

(^) The time of the maximum (or minimum) TMAX, and the maximum (or 
minimum) value ANSMAX is found by fitting a parabola to the three 
points T(I), ANS(I), for 1=1, 2, 3. 

(5) The threshold level TRLEV is computed as a fixed fraction of the 
maximum (or minimum) value ANSMAX 

I 

Input Parameters: 

INDEX1 Array index at which to start search for the maximum. 

T Array of time values 

ANS Array of electronics output signal values 

ISIGN Max or Min search option. If ISIGN is less than zero the 
program searches for a minimum. 

Threshold coefficient 


THRESH 


Output Parameters? 


IMAX 

TRLEV 

TMAX 

ANSMAX 


Array index nearest the maximum (or minimum) point 

Threshold level 

Time associated with maximum 

Estimate of maximum value of AMS array 


GUESS 


Purpose: The subroutine GUESS finds the initial estimate of the 

threshold crossing time for beginning the Newton-Raphson 
iterative procedure. 


Procedure; 

(1) The search procedure starts from array index associated with the 
maximum voltage point in the electronics output voltage array* 
ANS. 

(2) The array index is decreased until the signal value becomes less 
than threshold value (or greater than threshold value if ISIGN is 
less than zero). This array index provides the pointer to the 
initial time estimate for the Newton-Raphson iterative procedure. 

(3) Step 1 through 2 are repeated except that the array index 
increased. 

Input Parameters; 

TRLEV Threshold level 

ANS Array containing electronics output signal 

IMAX Array index of the maximum point 

ISIGN Positive or negative crossing search option. If ISIGN is 
less than zero the array indices are changed until the array 
value is above the threshold. 

Output Parameters; 

ICR0SS(2) Array indices nearby the threshold crossings 
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RTNI 


Purpose: The subroutine RTNI solves for the zeros of a function given 

an intial guess by using the Newton-Raphson iteration 
scheme. The external reference^ FCT provides the function 
value and its derivative to RTNI. 


FCT 

Purpose; FCT is a subroutine that calculates the value of the 
derivates whose zeros are found by the Newton-Raphson 
iterative procedure, RTNI. FCT evaluates the electronics 
output signal by calling AIR to interpolate on the tables of 
the output voltages vs. time. The function value returned 
, by FCT is the output signal minus the threshold level. The 

output signal derivative is evaluated numerically. 


SECTION 6 - SCANNER MEASUREMENT PREDICTOR 


6.1 FUNCTIONAL DESCRIPTION 

The Scannec Measurement Predictor (SMP) performs the following 
functions: 

(1) Simulates the sensor measurements from each Landsat-D conical 
scanner based on given ephemerls and reference attitude data. 

(2) Generates a dataset containing spacecraft position and attitude 
in^r .lation and one of the following 

Mode A. (Predicted/Observed). Predicted horizon scanner 
measurements along with actual measurements obtained 
from the Telemetry Data Base 

Mode B. (Predicted/Predicted). Measurements predicted from two 
sets of model parameters 

(3) Generates basic statistics on the information written to the 
output dataset. 

(i<) Generates a printout containing the predicted horizon crossing 
positions on the Earth. 

The output dataset is read by the Data Plotting and Fitting Utility in 
support of the sensor evaluation activities. 

To obtain the predicted measurements, the SMP first computes scanner 
measurements by modeling the intersection of the scan cone with an 
ellipsoid or oblate Earth assuming nominal conditions, and then 
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applies corrections bo these initial estimates, The correction 
factors simulate various error sources including misalignment of the 
scanner, responses of the sensor optics and electronics to variations 
in Earth horizon radiance, sensor output voltage and telemetry 
calibrations, sensitivity of the scanner measurements bo temperature 
changes, and presence of random noise. 

The two prediction modes A and B are for inflight data and theoretical 
analysis respectively. Mode A is the basic operational mode for 
inflight data analysis. It provides a dataset from which the residual 
errors between the observed and predicted measurements can be plotted 
and analyzed. Mode B is provided for convenient analysis of the 
effects of sensor model parameter variations. For example, this mode 
can be used to show the corrections to the scanner measurements due to 
effects of Earth oblabeneas and horizon radiance variations, More 
generally, this mode can be used to show the effects of any sensor 
model parameter changes that are considered for analysis purposes. 
After the output dataset is generated in this mode, the Data Plotting 
and Fitting Utility can be run to readily provide plots of the effects 
of the parameter changes as a function of time or other available 
parameter. Data fitting to the mode B output also provides an 
important system capability. In particular, the effects of systematic 
horizon radiance variation as predicted by the SOES utility and stored 
in the Triggering Heights Data Base can be fit by selected finite 
Fourier series parameters to provide an analytical representation of 
these effects. The systematic horizon radiance variation effects can 
be fit two ways, as a function of orbit position, and as a function of 
horizon crossing latitude. 

Figure 6-1 diagrams the basic procedure used in the SMP. Each of the 
steps in the diagram are described below. Reference 4 contains a full 
description of the mathematics of the model. 


6-2 


6-3 

















stop 1 
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The model parameters are read from a NAMELIST input file. The 
parameters Include data needed for simulation of error sources as well 
as control information that indicates which input/output options have 
been selected. Section 6.2.1 lists the NAMELIST parameters in detail. 

Step 2 


The following parameters, which do not change throughout the 
simulation, are initialized: 

1. CONE (1) “ scan cone radial arc length from scanner ^ spin axis 

2. AZIM (i) - scanner mounting alignment azimuth angle 

3. TILT (i) - soannerjL mounting alignment tilt angle 

A. TVIST(i) - scanner mounting alignment twist angle 

5. SHOUNT (3»3»i) - spacecraft body to scannerj^ coordinates 

transformation matrix 

The subscript 1 indicates that these parameters are needed for both 
scanner 1 and scanner 2, 

Step 3 

If the ephemeris or attitude data are obtained from the Telemetry 
Data Base, or If the output dataset is to contain predicted vn. 
observed measurements (as opposed to predicted vs. predicted 
measurements), the Telemetry Data Base is read in this step. 

One record is read from the Telemetry Data Base at a time, Each 
record contains 4 OBC reference attitudes, 4 ephemeris, and 128 
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scanner pitch and roll measurements. See Section 2.1.2 for a detailed 
description of the telemetry record. 

Step M 

The reference attitude and ephemerls data are Initialized, There are 
two possible sources for the reference attitudes: 

(1) Telemetry Data Base - These are the OBC attitudes that are 
Input as Euler parameters. 

(2) Constant Pitch, 5ioll, Yaw - The pitch, roll, and yaw 
parameters are obtained from the NAMELIST input file. 

There are also two sources for the ephemerls: 

(1) Telemetry Data Base - The telemetry contains the spacecraft 
position and velocity vectors in GCI coordinates. 

(2) Ground Ephemerls - There is a utility program called EPHEMX 
that reads a NASA standard ephemerls tape or uses NAMELIST 
Input orbital elements. 

If either the attitude or ephemerls data are obtained from the 
Telemetry Data Base, the system defaults to using the telemetry timing 
information also. In this case, there is an additional option to 
interpolate the attitude/ephemeris measurements, which are available 
at a rate of 4 samples per major frame, to enable a full comparison 
with the scanner measurements which arc sampled 128 times per major 
frame. If this option is selected, the interpolation is performed in 
this step by a simple linear method. 
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step 5 


The scanner measurements are Initially in the form of voltage outputs. 
The voltages are oonverted to counts when they are downlinked the 
telemetry stream. The Telemetry Prooessor stores the measurements in 
the Telemetry Data Base in this form. 

In this step, the soanner pitch and roll measurements are converted to 
angular Earth width and phase, and scanner and bolometer temperature 
measurements are oonverted to degrees centigrade. 

Step 6 


The following parameters must be initialized for every new measurement 
prediction: 


1. 

ATT(3,3) 

■ 

orbital to spacecraft body coordinates transforma- 
tion matrix (specified by pitch, roll, yaw) 

2. 

0RBNRMC3) 

- 

orbit normal direction 

3. 

GCIORB(3,3) 


geocentric- inertial to orbital coordinates 
transforlntition matrix 


GCISC(3,3) 


geocentric inertial to S/C Body coordinates 
transformation matrix 

5. 

ES(3,i) 

- 

scanner to Earth center direction 

6, 

ETA(i) 


arc-length separation between scanner^^ spin axis A 
and Earth vector E 

7. 

SUBLAT 

- 

latitude of subsatellite point 

8. 

RSUB 

- 

Earth radius at subsatellite point 

9. 

RHO(i) 

- 

Earth angular radius as seen from S/C 

10. 

SPINAX(3,i> 

- 

Scanner^ spin axis 
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The subscript i indicates that the parameter is dependent on the 
scanner number. 

Step 7 

A test is made to see whether or not the scan path crosses the Earth 
horizon. If not, the scanner measurements cannot be predicted. Flag 
values will bo entered in the output in this case. 

Stop 8 


The initial estimates of the Earth width and phase angles are 
calculated. The method used in the calculations involves iteratively 
computing horizon crossing positions on an ellipsoid model of the 
E.rth. The initial estimates of Earth width and phase are 

adjusted according to the following constant bias parameters; 

1. Constant adjustment to the horizon radiance height that triggers 
the sensor 

2. Constant adjustment to the Earth angular radius which models the 
effects due to finite sensor field-of-view 

3. Constant adjustments to the Earth width and phase angles which 
model effects due to sensor electronics (e.g., circuit response 
delays or errors in the amplifier gains). These corroctions are 
applied in the form of a polynomial with user-specified 
coefficients. 

Step 10 

In this step, the initial estimate for Earth width and phase are 
adjusted to account for sensor responses to systematic variations in 
Earth horizon radiance. The correction factors usually are obtained 
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from the Triggering Heights Data Base which is generated by the Sensor 
Optics and Electronics Simulator, An alternative to using the data 
base is to supply an analytic function that represents the contents of 
the data base. If no adjustments from the Triggering Heights Data 
Base are desired, this step can be skipped. 

Step 11 

Temperature dependent biases can be modelled by adding correction 
factors that are computed as a function of the temperature data 
available in the Telemetry Data Base, 

Step 12 


Steps 6 through 11 are performed twice, i,e, once for each scanner. 
If both scan paths do not cross the horizon,, steps 13 through 16 
described below should be skipped. 

Step 13 

If the comparison type is predicted vs. predicted, steps 6-12 must be 
repeated to obtain the second set of predicted measurements. 

Step in 

In the case in which the output dataset is to contain predicted 
measurements vs. predicted measurements, it is possible to simulate 
the presence of random noise in one set of measurements. This option 
enables the comparison of "ideal” measurements with simulated noisy 
measurements. 
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st ep 15 


The residual errors between the predicted measurements and observed 
measurements (or seoond set of predicted measurements) are calculated. 
All computations necessary to obtain the basic statistics for the 
summary report are made. The basic statistics include the mean 
residual error, root-mean-square residual error, and an estimate of 
the standard deviation of the residual errors as well as the minimum, 
maximum, and mean of the predicted and observed Earth width and phase 
measurements, scanner temperature, bolometer temperature, reference 
attitude, and spacecraft altitude. 

Step 16 

As each set of predicted measurements is computed. It is added to the 
output dataset. Note thot. steps 6 through 12 are performed twice, 
i.e. once for each scanner. Thus each output record contains all 
relevant measurements for both scanners. The SMP optionally generates 
a printout of the horizon positions viewed by the scanners as a 
'"unction of time. The Horizon Positions printout aids in the 
empirical evaluation of the effects of cold clouds on the scanner 
measurements. The computed horizon positions are based on the 
reference attitude and on the NAMELIST input parameters DHIN and 
DHOUT. The NAMELIST parameters IPRNT determines whether or not the 
Horizon Positons printout is generated. If IPRNT=1, a record of the 
printout is output in this jtep. 

Step 17 

If the attitude or ephemeris data are o’’tained from the Telemetry Data 
Base, there are either 4 or 128 time frames (a time frame is 
represented by one reference attitude and one ephemeris measurement) 
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extracted from each record of the telemetry dataset. Steps 6 through 
16 are repeated for each time frame. 

Step 18 

After all of the time frames generated from one telemetry record have 
been processed, steps 3 through 17 are repeated for each successive 
telemetry record. 


Step 19 

The final statistical computations are made and the summary report is 
generated. The summary report contains all of the model parameters, 
descriptive information about the simulation (e.g., time span 
covered), and the computed statistics. 
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6.2 


INPUTS AND OUTPUTS 


6.2.1 NAMELIST Inputs 

The NAMELIST input file for the Scanner Measurement Predictor, 
4SMPIN, contains the following variables: 

Control Parameters 


ICOHP Flag indicating type of comparison to be generated: 

ICOMP = 0 Observed vs. predicted 
= 1 Predicted vs, predicted 

= 2 Predicted vs. predicted with noise simulated 
in second set of measurements 
IFIRST First record of Telemetry input file to be read 
LAST Last record of Telemetry input file to be read 

TSTART Start time to read from Telemetry input file in format 

YYMMDD.HHMMSS 

TEND End time to read from Telemetry input file in format 
YYMMDD.HHMMSS 

TSTEP Time step for predicted vs. predicted mode comparisons 

NP Flag indicating number of predictions to be made when 

using timing information from the telemetry dataset: 

= 1 make 1 prediction per telemetry record. 

= make prediction s per telemetry record. 

= 128 make 128 predictions per telemetry record. 

Reference Attitude Parameters 

lATT Flag for source of reference attitudes: 

lATT = 0 OBC Attitudes from Telemetry Data Base 
= 1 Constant Pitch, Roll, Yaw 
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PITCH Pitch angle for constant attitude 
ROLL Roll angle for constant attitude 
YAW Yaw angle for constant attitude 

Ephemerls Parameters 

lEPHEM Ephemeris source flag; 

lEPHEM = 0 Use Telemetry Data Base 

= 1 Use ephemeris subroutine EPHEMX 

ISPC EPHEMX methods for obtaining spacecraft positon. 

= 1 OREHEN, use orbital elements from NAMELIST 
r, 2 DTAPRE (entry point of ROITAP), read 
sequential EPHEM file 

= 3 GETVCT/GETHDR read direct access orbit file 


TEPOCH 

Epoch time of orbital elements in 

format YYMMDD.HHMMSS 

A 

Semi-major axis (km) 


E 

Eccentricity 


EYE 

Inclination (degrees) 


EMO 

Mean anomaly (degrees) 


WO 

Argument of perigee (degrees) 


RANODE 

Right ascension of ascending node 

(degrees) 


Scanner Alignment Parameters 

AZIM0(2) Scanner mounting alignment azimuth angle 
TILT0(2) Scanner mounting alignment tilt angle 
TWIST0(2) Scanner mounting alignment twist angle 
DAZIMC2) Misalignment in azimuth angle 
DTILT(2) Misalignment in tilt angle 
DTWIST(2) Misalignment in twist angle 
CONEO(2) Nominal scanner cone angle 
DC0NE(2) Cone angle bias 
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Earth ModellnR Parameters 


ERAD Earth equatorial radius 

FLATCO Earth flattening coefficient 

DHIN Degrees of scanner rotation Inside the Ideal Earth-ln 
crossing used In the horizon crossing Earth latitude and 
longitude determination. 

DHOUT Degrees of scanner rotation inside the Ideal Earth-out 
crossing used In the horizon crossing Earth latitude and 
longitude determination. 

Horizon Radiance Correction P arameter s 


IHRCOR Flag for horizon radiance corrections; 

IHRCOR s 0 Use Horizon Triggering Heights Data Base 
= 1 Use analytic expression 
= 2 Do not apply corrections 

HORHGT Nominal horizon radiance height for input to oblate 
Earth model 

DHRHGT Constant horizon radiance height bias 
IVAR Flag indicating independent variable in analytic 
expression; 

IVAR = 0 horizon crossing latitude 

= 1 orbit position (true anomaly from the 
ascending node) 

MXFOUR Maximum order of Fourier expansion (maximum of 7) 

H Zero-order coefficient in Fourier expansion 

A Sin-wave coefficients 

B Cos-wave coefficients 

HRNOM Nominal horizon radiance height for correction factor 


Temperature Dependent Corrections 


ITPCOR 


TNOM(2) 

BTN0M(2) 

ETCOEF 

(3,2) 

HTCOEF 

(3,2) 

EBCOEF 

(3,2) 

HBCOEF 

(3,2) 


Flag for •> ■'mperature dependence corrections: 

ITPCOR = 0 Do not apply corrections 
= 1 Apply corrections 
Nominal scanner temperature 
Nominal bolometer temperature 

Coefficients in polynomial correction for scanner 
temperature dependence in E-voltage output 
Coefficients in polynomial correction for scanner 
temperature dependence in H-voltage output 
Coefficients in polynomial correction for bolometer 
temperature dependence in E-voltage output 
Coefficients in polynomial correction for bolometer 
temperature dependence in H-voltage output 


Constant Bias Parameters 


DERAD(2) Constant angular Earth radius bias 
WIDTH0(2) Nominal Earth width 
PHASE0(2) Nominal Earth phase 

WC0EFC2) Coefficients for constant bias and calibration 
polynomial correction to Earth width prediction 
PC0EF(2) Coefficients for constant bias and calibration 
polynomial correction to Earth phase prediction 

Voltage Conversion Parameters 

C0NWID(2) Conversion factor from Earth width to rotation angle Eij. 
C0NPHA(2) Conversion factor from Earth phase to rotation angle 
EVT0A(2) E voltage-to-angle scale factor 
HVT0A(2) H voltage-to-angle scale factor 


6-14 



DPITCHC2) 
DROLL (2) 
CTOV 
VTOST 
VTOBT 


Pitch offset voltage 

Roll offset voltage 

Counts-to-volts scale factor 

Voltage to scanner temperature conversion factor 

Voltage to bolometer temperature conversion factor 


Noise Simulation Parameters 


NPOINT Number of points used in N-polnt average for noise 
simulation 

STDE Standard deviation of noise in E voltage 

STDH Standard deviation of noise in H voltage 

Input/Output Control Parameters 

ITLM Logical unit number of telemetry input file 

ITHDB Logical unit number of Horizon Triggering Heights input 

file 

IDBG Logical unit number of debug output file 

IMDB Logical unit number of predictor output dataset 

IHPP Logical unit number of Horizon Positions Printout 

IPRNT Flag indicating whether or not to generate the Horizon 
Position Printout; 

IPRNT =0 do not generate printout 
= 1 generate printout 

XHPSFT Skipping factor Tor Horizon Positions Printout (output 
every n^*’ time frame) 

LEVDBG Level of debug output to be generated: 

LEVDBG = 0 none 

= 10 minimal 
= 20 moderate 
= 30 maximum 
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Dataset Inputs 


Telemetry Input 


The telemetry dataset is selected from t.he Telemetry Data Base, See 
Section 2.1 for a description of its format. 

TriKKering Heights Input 

The Triggering Heights Data Base is read if the option to apply 
horizon radiance corrections from it is selected. See Section 2 .^ for 
a description of its content and format, 

6.2.3 Dataset Outputs 

The Scanner Measurement Predictor writes to a dataset that is entered 
In the Measurements Data Base. See Section 2.5 for a description of 
its content and format, 

6.2.4 Printouts 
Summary Report 

f * 

The summary report contains the following information: 

(1) NAMELIST Parameters 

(2) Time span covered in the simulation 

(3) Total number of measurements predicted 

(4) Minimum, maximum, mean, root-mean-square, standard deviation of 
the residuals 


(5) Maximum, minimum, and mean of eacb of the following variables: 

- Earth width and phase (predicted and observed) 

- Earth-in and Earth-out angles (predicted and observed) 

- Bolometer and scanner temperatures 

Horizon Positions Printout 


The Horizon Positions Printout displays the computed horizon crossing 
positons for each scanner as a function of time and orbit position. 
The horizon crossing positions (in-crossing and out-crossing) are 
expressed in terms of latitude and longitude. The orbit position is 
expressed in terms of subsatellite latitude and longitude. The time 
is in format YYMMDD.HHMHS3. 

Debug Output File 

The amount of information contained in the debug output file is 
determined by the NAMELIST parameter LEVDBG. The actual output 
associated with each level has not been defined. 


i 


!» 

I 

f 

f- : 

6.3 SYSTEM STRUCTURE f 

Figure 6~2 shows the baseline diagram of the Scanner Measurement ^ 

Predictor. The MAIN routine acts as the driver of the system and 
calls the remaining subroutines to perform moat of the computation. 

Table 6-1 lists the functional steps described in Section 6.1 along 
with the names of the major subroutines that support them. 
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Table 6-1 Relationship Between Funotlonal Steps and Subroutines 


4 


Step 

Subroutines 

1 

NAMLST 

2 

INIT 

3 

RDTLM 

4 

ATTEPH 

5 

CONVRT 

6 

SETPRM 

7 

CHKPTH 

8 

PREDCT 

9 

BCORR 

10 

HRCORR 

11 

TCORR 

12 

MAIN 

13 

MAIN 

14 

NOISE 

15 

STATS 

16 

MAIN 

17 

MAIN 

18 

MAIN 

19 

STATS, SUHHRY 
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6.4 


SUBROUTINE DESCRIPTIONS 


Eaoh subroutine of the Scanner Measurement Predictor subsystem is 
described below in terms of its purpose, procedure, input parameters, 
and output parameters. Descriptions of existing software (o.g. 
TC0N40, UNVEC) are omi: ted. 

MAIN 


Purpose: MAIN is the driver of the Scanner Measurement Predictor. 

Procedure; 

(1) Call NAMLST to input model parameters. 

(2) Call INIT to initialize orbit-level parameters that do nob 
change throughout the simulation, 

<3) If lATTaO or lEPHEMaO or IC0HP=0, call RDTLM to read a 
record of the telemetry dataset. 

(4) Call ATTEPH to initialize attitude and ephemeris dependent 
parameters. 

<5) If IC0MP=0 or ITPCORal, call CONVRT to convert scanner 
measurements from counts to angles and °C. 

(6) Call SETPRM to initialize scanner dependent parameters. 

(7) Call CHKPTH to determine whether or not the scan path 
crosses the ho^^ ion. If not, insert a flag for the 
predicted measurements and omit steps 8 through 11. 

(8) Call PREDCT to calculate Earth width and phase, WIDTHP and 
PHASE P. 

(9) Call BCORR to apply constant bias corrections to WIDTHP and 
PHASER. 
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(10) If IHRC0R=1, call HRCORR to apply horizon radianoe variance 
corrections to WIDTH? and PHASE?, 

(11) If ITPC0R=1, call TCORR to apply temperature dependent 
corrections to WIDTHP and PHASEP. 

(12) Loop through steps 6 to 11 for each scanner, 

(13) If ICOHPsO, repeat steps 6 to 12 to obtain second set of 
predicted measurements. 

(14) If IC0HP=2, call NOISE to simulate noise in second set of 
predicted measurements. 

(15) Call STATS to compute residuals and to update statistical 
computations. 

(16) Write record to output dataset (Measurements Data Base). 
See Section 2.5,2 for a description of the record contents. 
If Horizon Positions Printout is requested, write out a 
record. 

(17) Afte'' all requested data has been processed, generate 
summary report. 

(18) Throughout the above steps, MAIN performs error checking and 
writes debug output messages if requested. 

Input Parameters; None 

Output Parameters: None 
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NAMLST 


Purpose: NAMLST reads in the model parameters from the NAMtiLIST input 

file &SMPIN. 


Procedure: 

O) Perform NAMELIST read. 

(2) Test for read errors and set error flags if any occurs. 

Input Parameters: None 

Output ParaiO''ters; 

• lERR Error flag; 

lERRsO no error 

=1 read error occurred 

The remaining output parameters are listed by common block: 

/CTRL/ 

ICOMP Comparison type flag: 

i 

ICOMP = 0 Predicted vs. Observed 

f 

= 1 Predicted vs. Predicted 

= 2 Predicted vs. Predicted (Noisy) 

NP Flag indicating number of predictions per major frame: 

NP = 1 1 prediction per major frame 

= 128 128 predictions per major frame 

=4 4 predictions per major frame 

TSTEP Time step to be used when telemetry timing data is not 
used. 
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/TLMCTL/ 


IFIRST 

LAST 

TSTART 

TEND 


/ATTIT/ 

lEPHEM 


lATT 


PITCH 

ROLL 

YAW 

/PARAMS/ 

AZIMO 

TILTO 

TWISTO 

DAZIM 

DTILT 

DTWIST 

CONEO 

DCONE 

ERAD 


First record of Telemetry input file to be read 

Last record of Telemetry input file to be read 

Start time to read from Telemetry input file in format 

YYMHDD.HHMMSS 

End time to read from Telemetry input file in format 
YYMHDD.HHMMSS 


Flag indicating source of ephemerls: 
lEPHEM = 0 Telemetry data base 

= 1 EPHEMX ephemeris model 
Flag for source of reference attitudes: 
lATT = 0 OBC Attitudes from telemetry data base 
= 1 Constant Pitch, Roll, Yaw 
Pitch angle for constant attitude 
Roll angle for constant attitude 
Yaw angle for constant attitude 


Scanner mounting alignment azimuth angle 
Scanner mounting alignment tilt angle 
Scanner mounting alignment rw^st angle 
Misalignment in azimuth angle 
Misalignment in tilt angle 
Misalignment in twist angle 
Normal scanner cone angle 
Cone angle bias 
Earth equatorial radius 
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FLATi.0 

HORHGT 

DHIN 

DHOUT 

/HRCORR/ 

IHRCOR 


HRNOH 

MXFOUR 

H 

A 

B 

DHRHGT 

IVAR 

/TCORR/ 

ITPCOR 


TNOM 

BTNOM 

ETCOFF 


Earth flattening coeffioient 

Nominal horizon radiance height (input to oblate Earth 
model) 

Earth-in horizon crossing position determination 
adjustment . 

Earth-out horizon crossing position determination 
adjustment. 


Flag for horizon radiance corrections: 

IHRCOR = 0 Use Horizon Triggering Heights Data Base 
= 1 Use analytic expression 

=2 Do not apply corrections 

Nominal horizon radiance height 
Maximu'^ order of Fourier expansion 
Zero-order coefficient in Fourier expansion 
Sin-wave coefficients 
Cos-wave coefficients 
Constant horizon radiance height bias 
Independent variable selection flag 


Flag for temperature dependence corrections: 

ITPCOR =0 Do not apply corrections 

- 1 Use scanner and bolometer temperatures for 
corrections 

Nominal scanner temperature 
Nominal bolometer temperature 

Coefficients in polynomial correction for scanner 
temperature dependence in E-voltage output 
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HTCOEF Coefficients in polynomial correction for scanner 
temperature dependence in H-voltage output 
EBCOEF Coefficients in polynomial correction for bolometer 
temperature dependence in E-voltage output 
HBCOEF Coefficients in polynominal correction for bolometer 
temperature dependence in H-voltage output 

/BIASES/ 

DERAD Constant angular Earth radius bias 
WIDTHO Nominal Earth width 
PHASED Nominal Earth phase 

WCOEF Coefficients for polynomial correction to Earth width 

PCOEF Coefficients for polynomial correction to Earth phase 

/CONVRT/ 

CONWID Conversion factor from Earth width to rotation angle 

CONPHA Conversion factor from Earth phase to rotation angle H.p 

EVTOA E voltage-to-angle scale factor 

HVTOA H voltage-to-angle scale factor 

DPITCH Pitch offset voltage 

DROLL Roll offset voltage 

CTOV Counts-to-volts scale factor 

/NOISE/ 

NPOINT Number of points used in N-point average for noise 
simulation 

SIDE Standard deviation of noise in E voltage 

STDH Standard deviation of noise in H voltage 
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/LUNITS/ 


ITLM Logical unit number of telemetry input file 
THDB Logical unit number of Horizon Triggering Heights input 
file 

IMDB Logical unit number of predictor output dataset 

IHPP Logical unit number of Horizon Positions Printout 

IPRNT Horizon Positions Printout on/off flag 

/DEBUG/ 

IDBG Logical unit number of debug output file 
LEVDBG Level of debug output to be generated 


original paCj. 

OF POOR QUAUTV 


INIT 


Purpose: INIT initializes orbit-dependent parameters that do not 

change throughout the simulation. 


Procedure; 

(1) Compute scanner mounting alignment angles; 
AZIHCI) = AZIMO(I) + DAZIHd) 

TILT (I) = TILTO(I) + DTILT(I) 

TWIST(I) = TWISTO(I) + DTWIST(I) 
CONE(I) = CONEO(I) + DCONE(I) 
for I = 1, 2 


( 2 ) 


Compute scanner mounting alignment matrix; 


SHOUNTd.I.I) 

SM0UNT(2,1,I) 

SMOUNT(3,1,I) 

SM0UNT(1,2,I) 

SM0UNT(2,2,I) 

SMOUNT(3,2,I) 

SM0UNT(1,3,I) 

SMOUNT(2,3,I) 

SMOUNT(3,3,I) 


s oos<TILT(I))»cos(AZIM(I)) 

= -cos(TWIST(I))*3in(AZIM(I)) 

+ (TILT(I))«cos(AZIM(D) 

= sin(TWIST(I))«sin(AZIM(D) 

+ cos(TWIST(I))«siri(TWIST(I))»cos(AZIM(I)) 
= cos(TWIST(I))»sin(AZIM(D) 

= cos(TWIST(I))»cos(AZIM(I)) 

+ 3in(TWIST(I))*sin{TILT(I))»sin(AZIM(D) 
= -sin(TWIST(I))»cos{A4lMa)) 

+ cos(TWIST{I))*3in(TILT(I))»sin(AZIH(I)) 
= -sin(TILT(D) 

= sin(TWIST(I))*cos(TILT(D) 

= co3(TWIST(I))«cos(TILT(I)) 


for I = 1,2 
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Input Parameters: 


ICOMP Comparison type flag 

The remaining i^^put parameters are from cpmmon block /PAPAMS/. 
Output Parameters: 

AZIM Scanner mounting alignment azimuth angle 

TILT Scanner mounting alignment tilt angle 

TWIST Scanner mounting alignment twist angle 

CONE Scan cone radial arc-length from scanner spin axis 

SMOUNT Scanner mounting alignment matrix 
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RDTLM 


Purpose; RDTLM reads a record of the telemetry input dataset. If 
IFIRST or TSTART is set, RDTLM positions the dataset at the 
appropriate record by performifis sequential reads. 


Procedure: 

(1) If IFIRST is greater than 1, perform sequential reads to 
position the dataset at record number IFIRST. 

(2) If ISTART is greater than 0, perform sequential reads until 
the time value from the record is greater than or equal to 
TSTART, 

(3) Read in the telemetry record. 

(4) Check for end-of-file conditions; 

a) Physical EOF 

b) Current record greater than LAST 

c) Time on current record greater than TLAST 

(5) Chock for read errors. 

Input Parameters; 

ITLM Logical unit number of the telemetry dataset 
IREC Number of the record currently being processed 
TCl'RR Time value of the rec^ord currently being processed 

The following input parameters are from corrjiion block /TLMCTL/: 

IFIRST First record of Telemetry input file to be read 
LAST Last record of Telemetry input file to be read 
TSTART Start time to read from Telemetry Input file in format 
YYHMDD.HHMMSS 
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1 / 


TEND End time to read from Telemetry input file in format 
YYMMSS.HHMMSS 

Output Parameters; 

lERR End-of-file and error flag: 
lERR = 0 successful read 
= 1 end-of-file 
s 2 read error 

IREC, TCURR are updated if successful read occurred 

The remaining output parameters are listed by common block; 

/CSDATA/ 

IROLLF Observed i'' 11 measurements in counts 
IPITF Observed pitch measurements in counts 
ISIGNL Signal status flags 

ITEMP Observed scanner temperatures in counts 
IBOLO Observed bolometer temperatures in counts 
ISENSR Sensor status flags 

MODES Mode flags indicating various operational conditions 
/TIMING/ 

TIME Time on each record in seconds since Sept, 1, 1957. 
/ATTEPH/ 

EPARMS Euler parameters representing OBC reference attitudes 
EPHPOS Spacecraft position in GCI coordinates 
EPHVEL Spacecraft velocity in GCI coordinates 
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ATTEPH 


Purpose: ATTEPH obtains the attitude and ephemeris data for a set of 

time frames. The number of time frames that are Initialized 
on eaoh call to ATTEPH depends on the NAMELIST parameters 
ICOHP and NP. 


Procedure: 

(1) Call GETATT to obtain the reference attitudes. If NP=128, 
and IATT=0, call INTERP to interpolate the measurements 
returned from GETATT. 

(2) Call GETEPH to obtain the ephemeris. If NP = 128 and 
IEPHEM=0, call INTERP to interpolate the measurements. 

(3) Perform error checking. 

Input Parameters: 

lATT Flag indicating attitude source 

lEPHEM Flag indicating ephemeris source 

NP Flag indicating number of time frames initialized 

TCURR Current time 

The remaining input parameters are listed by common block: 

/ATTEPH/ 

EPAiilS Euler parameters representing OBC reference attitudes 
EPHPOS Spacecraft position in GCI coordinates 
EPHVEL Spacecraft velocity in GCI coordinates 
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/DEBUG/ 


IDBC Logical unit number of debug output file 
LEVDBG Level of debug output to be generated 

Output Parameters: 

XPITCH Pitch measurements 

XROLL Roll measurements 

XYaW Yaw measurements 

SCPOS Spacecraft position vectors 

SCVEL Spacecraft velocity vectors 

lERR Error flag: 

lERR = 0 no error 

= 1 error occurred 
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GETATT 


Purpose: GETATT obtains reference attitude in the form of pitch* 

roll, yaw from the OBC Euler parameters. 

Procedures: 

(1) Compute the attitude matrix from the Euler parameters. 

(2) Obtain pitch, roll, and yaw from the attitude matrix. 

Input Parameters; 

EPARHS Euler parameters reprersnttng OBC reference attitudes 

Output Parameters: 

PITCH Pitch measurements 
ROLL Roll measurements 

YAW Yaw measurements 
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GETEPH 


Purpose: GETEPH obtains the ephemeris in the form of spacecraft 

position and velocity vectors in GCI coordinates. 

! 

Procedures: ! 

I I 

i 

(1) If lEPHEH = If call EPHEHX to obtain position and velocity 
vectors, 

(2) If lEPHEM = 0, use ephemeris from telemetry. j 

j 

' ( 

Input Parameters: ,| 

I [ 

Ephemeris source flag 

Number of ephemeris points to be retrieved (1 or 4) 

Spacecraft velocity vector from Telemetry Data Base 
Spacecraft position vector from Telemetry Data Base 
Current time 
Time increment 

Output Parameters: 

E Spacecraft to Earth center vector 

V Spacecraft velocity vector 





lEPHEM 

NP 

EPHVEL 

EPHPOS 

CURTIM 

TIMSTP 
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CONVRT 


Purpose; CONVRT converts the observed pitch and roll measurements 
from telemetry counts to angular Earth width and phase, and 
converts the observed temperatures from telemetry counts to 
degrees centigrade. 


Procedures: 

(1) Convert voltages to angular Earth width and phase: 
WIDTH(I,J)=1/C0NWID»(1/EVT0A*(VPITCH(I,J)-DPITCH)) + WIDTHO 
PHASEd, J)=1/C0NPHA*(1/HVT0A*(VR0LL(I,J)-DR0LL)) + PITCH 
for I = 1,128; J = 1,2. 

(2) Compute Earth-in and Earth-out angles; 

AIN = WIDTH/2 + PHASE 

AOUT = -WIDTH/2 + PHASE 

(3) Convert scanner and bolometer temperatures: 

TEMP(I) = ITEMPCI, J)/CT0V » TVTOD 

BOLD (I) = IB0L0(I,J)/CT0V * BVTOD 
for I = 1,128; J = 1,2. 

Input Parameters: 

IFLAG Flag indicating conversions to be pe'-fr»"med: 

IFLAG = 0 all 

= 1 temperature only 

The remaining input parameters are listed by common blocks: 


/CONVRT/ 


CONWID Conversion factor from Earth width to rotation angle E-p 

CONPHA Conversion factor from Earth phase to rotation angle H-j- 

EVTOA E voltage-to-angle scale factor 

HVTOA H volbage-to-angle scale factor 

DPITCH Pitch offset voltage 

DROLL Roll offset voltage 

CTOV Counts-to-volts scale factor 

/CSDATA/ 

IROLLF Observed roll measurements in counts 
IPITF Observed pitch measurements in counts 
ISIGNL Signal status flags 

ITEMP Observed soanre* temperature in counts 

IBOLO Observed temperature in counts 

ISENSR Sensor status flags 

MODES Model flags indicating various operational conditions 
Output Parameters: 


WIDTH 

Observed 

Earth width measurements 

PHASE 

Observed 

Earth phase measurements 

AIN 

Observed 

Earth~in angle 

AOUT 

Observed 

Earth-Out angle 

TEMP 

Observed 

scanner temperature (deg C) 

aoLO 

Observed 

bolometer temperature (deg C) 
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SETPRM 


Purpose: SETPRM Initializes parameters that are dependent on each 

bime-frame and scanner. 

Procedures: 

(1) Compute orbit normal direction. 

(2) Compute GCI to orbital coordinates transformation matrix. 

(3) Compute GCI to S/C coordinates transformation matrix. 

Compute scanner to Earth center direction. 

(5) Compute arc-length separation between scanner spin axis and 

Earth vector. 

(6) Compute latitude of subsatellite point, 

(7) Compute Earth radius at subsatellite point. 

(8) Compute Earth angular radius as seen from S/C. 

(9) Compute scanner spin axis. 

Input Parameters: 

V Spacecraft velocity vector 

E Spacecraft to Earth center vector 

ATT Attitude matrix 

SMOUNT Scanner mounting alignment matrix 
ERAD Earth equatorial radius 
FLATCO Earth flattening coefficient 
HORHGT Nominal horizon radiance height 

Output Parameters: 

ORBNRM Orbit normal direction 

GCIORB GCI to orbital coordinates transformaion matrix 

GCISC GCI to spacecraft body coordinates transformation matrix 
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■J 

ES Scanner to Earth center direction 

ETA Arc-length separation between scanner spin axis and Earth 
vector 

SUBLAT Subsatellite latitude 
RSUB Earth radius at subsatellite point 
RHO Earth angular radius 
SPINAX Scanner spin axis 
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CHKPTH 


Purpose: CHKPTH determines whether or not the scan path crosses the 

Earth horizon. 


Procedure; 

(1) Determine whether scan path intersects the Earth for at 
least part of the scan: 

ABSCETA - CONE) < RHO 

(2) If condition (1) is met, determine whether scan path crosses 
the horizon: 

(a) If ETA + CONE ^180, then 

ETA < CONE + RHO 

(b) If ETA + CONE >180, then 

ETA + CONE < 360 - RHO 

(3) If conditions (1) and (2) are met, set IPATH = 1. Otherwise 
set IPATH = 0. 

Input Parameters; 

ETA Arc-length separation between scanner spin axis and Earth 
vector 

CONE Scan cone radial arc-length from scanner spin axis 
RHO Earth angular radius as seen from the S/C. 

Output Parameters: 


IPATH Flag indicating whether or not scan path crosses horizon: 
IPATH = 0 path does not cross 
= 1 path crosses 
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PREDCT 


Purpose; PREDCT computes the predicted Earth width and phase 
measurements assuming an oblate model of the Earth. 


Procedure; 

(1) Compute initial approximation to Earth width and phase; 
WOs2<‘ACOS{(COS(RHO)-COS(ETA)COS(CONE))/(SIN(ETA)»SIN(CONE))) 
P0=ATAN2(-ES(2),ES(3)) 

(2) Call RFALSl to iteratively compute Earth-ln angle AIN with 
AINO = WO/2 + PO, and AIN1 based on sign and value of AINO, 

(3) Call RFALSl to iteratively compute the Earth-out angle, AOUT with 
AOUTO = -W/2 +P0, and AOUT1 based on sign and value of AOUTO. 

(i<) Compute final predictions of Earth width and phase by 
WIDTHP = AOUT - AIN 
PHASEP = (AIN + AOUT)/2 

Input Parameters: 

RHO Earth angular radius 

ETA Arc-length separation between scanner spin axis and Earth 
vector , 

CONE Scan cone radial arc-length 

ES Scanner to Earth vector 

E Earth vector 

ERAD Earth equatorial radius 

SMOUNT Scanner mounting alignment matrix 

GCISC GCI to spacecraft coordinates transformation matrix 

EPS Tolerance for convergence interval 

HORHGT Nominal horizon radiance height 
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Output Parameters; 


WIDTHP Earth width estimate 
PHASEP Earth phase estimate 
AIN Earth-in angle 
AOUT Earth-out angle 
lERR Error flag: 

lERR = 0 no error 

= 1 no convergence 


BCOAR 


Purpose: BCORR applies constant bias parameters to the initial Earth 

width and phase measurements. 


Procedure: 

(1) Apply constant angular Earth radius adjustment; 

WIDTH? s WIDTH? + (DWDR « DERAD) 

PHASEP = PHASE? + (DPDR » DLRAD) 

Where DWDR and DPDR are the partial derivatives of width and 

phase with respect to the angular Earth radius, 

(2) Apply polynomial correction; 

WIDTH? 5 WIDTH? + (WCOEF(l) + WC0EF(2) » (WIDTHP - WIDTHO) 

+ WC0EFC3) * (WIDTHP - WIDTHO) »«2 
+ WCOEF(4) » (PHASEP - PHASED) 

+ WC0EF(5) * (PHASFP - PHASED) »*2) 

PHASEP = PHASEP + (PCOEF(I) + PCOEF(2) * (PHASEP - PHASED) 

+ PCOEFO) " (PHASEP - PHASED) «»2 
+ PCOEF(4) » (WIDTHP - WIDTHO) 

+ PCOEF(5) * WIDTHP - WIDTHO) **2 


Input Parmeters: 

WIDTHP Earth width estimate 

PHASEP Earth phase estimiate 

DERAD Earth radius constant bias 

WIDTHO Nominal Earth width 
PHASED Nominal Earth phase 

WCOEF Earth width polynomial correction coefficients 

PCOEF Earth phase polynomial correction coefficients 
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Output Parameters: 


WIDTHP Earth width estimate 
PHASEP Earth phase estimate 
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RFALSI 


Purpose; 


RFALSI iteratively computes the solution of a nonlinear 
equation using the Regula Falsi method. 


Procedure; 

(1) Compute {n+1)st point by 

XNPl = XN - FN * ((XN-XNM1)/(FN-FNM^)) 
where FN * FNM1 < Or 

(2) Perform convergence test; 

ABS{FNP1) < EPS 

(3) If no convergence, iterate with 
XNM1 = XN 

XN = XHP1 
FNM1 = FN 
FN = FNPl 


Input Parameters; 

XO.X1 Initial points 

F External function for calculation of the equation at each 

X. 

EPS Convergence tolerance 

Output Parameters; 

X Root of the equation within the given tolerance 

lERR Error flag; 

lERR = 0 no error 

= 1 no convergence 
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HnCORR 


Purpose; HRCORR applies corrections for scanner responses to 
variations in Earth horizon radiance to the Earth width and 
phase estimates. 


Procedure; 

(1) If IHRCOR = 0, call RDTHD3 to obtain the correction factor 
DHT. 

(2) If IHRCOR s 1, evaluate analytic expression to obtain 
correction factor DHT 

(3) Apply constant bias correction factor DHRHGT. 

Input Parameters; 

IHRCOR Horizon radiance correction flag 
DHRHGT Constant bias horizon radiance height 
MXFOUR Order of Fourier expansion 
H Zero-th order term 

A Sin-wave coefficients 

B Cos-wave coefficients 

HRNOM Nominal horizon radiance height 
WIDTHP Earth width estimate 
PHASEP Earth pliase estimate 
IVAR Independent variable selection flag 

Output Parameters; 

WIDTHP Earth width estimate 
PHASEP Earth phase estimate 
lERR Error flag 

lERR = 0 no error 

= 1 error in reading THDB 
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TCORR 


Purpose; TCORR adjusts the Earth width and phase estimates to account 
for temperature dependent effects. 


Procedure: 

(1) Compute scanner temperature correction factor for Earth width; 
WIDTHP = WIDTHP + [ETCOEF(I) + ETC0EF(2) * (TEMP - TNOH) 

+ ETCOEF(3) * (TEMP - TNOM) . 

(2) Compute scanner temperature correction factor for Earth phase; 
PHASEP = PHASEP + [HTCOEF(I) + HTC0EF(2) “ (TEMP - TNOM) 

+ HTOCEF(3) » (TEMP - TNOM) »«2] 

(3) Compute bolometer temperature correction for Earth width: 
WIDTHP = WIDTHP + [EBCOEF(D + EDC0EF(2) * (BOLO - BTNOH) 

f: EBCOEF(3*) » (BOLO -* BTNOM) »*2] 

'4) Compute bolometer temperature correction for Earth phase: 
PHASEP = PHASEP + [HBCOEF(I) + HBC0EF(2) * (BOLO - BTNOM) 

+ HBCOEF(3) * (BOLO - BTNOM) »»2] 

Input Parameters; 

WIDTHP Earth Width estimate 

PHASEF Earth Phase estimate 

* 

TEMP Observed Scanner temperature 

TNOM Nominal Scanner temperature 

BOLO Observed bolometer temperature 

BTNOM Nominal bolometer temperature 

ETCOEF Coefficients for scanner temperature correction 
to Earth Width 

HTCOEF Coeificients for scanner temperature correction 
to Earth Phase 
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EDCOEF Coeffioients for bolometer temperature) correction 
to Earth Width 

HBCOEF Coefficients for bolometer temperature correction 
to Earth Phase 

Output Parameters; 

WIDTHP Earth Width estimate 

PHASEP Earth Phase estimate 


STATS 


Purpose: Compute basic statistics (minimum, maximum, mean, root-mean- 

square, standard deviation) of Earth width and phase 
measurements. Earth-in and Earth-out angles, residuals, 
temperatures. 

Procedures: 

(1) If IFLAG = 0, then compute residuals, update running sums, 
and update maxima and minima. 

(2) If IFLAG = 1, then compute final statistics. 

Input Parameters: 

IFLAG Flag indicating termination of data. 

IFLAG = 0 update Eums, c^'c, 

= 1 compute final statistics 
WIDTHP Earth Width estimate 

PHASEP Earth Phase estimate 

AINP Earth-in angle estimate 

AOUTP Earth-out angle estimate 

WIDTH Observed Earth Width 

PHASE Observed Earth Phase 

AIN Observed Earth-in Angle 

AOUT Observed Earth-out Angle 

TEMP Observed Scanner uempenture 

BOLO Observed Bolometer temperature 
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Output Parameters 


XMIN 

XMAX 

XMEAN 

RMS 

STDR 

RMIN 

RMAX 

RMEAN 

N 


Minima of above variables 

Maxima of above variables 

Means of above variables 

Root-mean-squares of residual errors 

Standard deviations of residual errors 

Minima of residual errors 

Maxima of residual errors 

Means of residual errors 

Total number of points 
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NOISE 


Purpose: NOISE adjusts the computed Earth width and phase estimates 

to simulate random noise. 

Procedure; 


(1) 

Call 

RAVNZ1 

to 

compute noise 

factor 

for 

Earth width estimate 


for 

scanner 

1. 






(2) 

Call 

RAVNZ2 

to 

compute noise 

factor 

for 

Earth width estimate 


for 

scanner 

2, 






<3) 

Call 

RAVNZ3 

to 

compute noise 

factor 

for 

Earth phase 

estimate 


for 

scanner 

1. 






(4) 

Call 

RAVNZ4 

to 

compute noise 

factor 

for 

Earth phase 

estimate 


for 

scanner 

2. 







Input Parameters: 

WIDTHP Earth Width estimate 

PHASEP Earth Phase estimate 

SIDE Standard deviation ol' noise in E voltage 

STDH Standard deviation of noise in H voltage 

NPOINT Number of points in running average 

Output Parameters: 

WIDTHP Earth Width estimate 

PHASEP Earth Phase estimate 


SUMMRY 


Purpose: SUMMRY generates a report of the prediction including basic 

statistics of various data. 


Procedure: 

Output each of the following; 

(1) Heading and date 

(2) NAMELIST parameters 

(3) Initial/ final times 

(4) Total number of measurements 

(5) Basic statistics 

Input Parameters; 

TINIT Initial time 

TFINAL Final time 

XHIN Minima of observed and computed measurements 

XMAX Maxima of observed and computed measurements 

XMEAN Means of observed and computed measurements 
RMS Root-mean-squares of residual errors 

STDR Standard deviations of residual errors 

RMIN Minima of residual errors 

RMAX Maxima of residual errors 

RMEAN Means of residual errors 

N Total number of points 

The remaining input parameters are from common blocks CTRL, TLMCTL, 
ATTIT, PARAMS, HRCORR, TCORR, BIASES, CONVRT, NOISE, LUNITS, DEBUG. 
See subroutine NAMLST for common block contents. 

Output Parameters; NONE 


/j 
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SECTION 7 - DATA PLOTTING AND FITTING UTILITY 


7.1 FUNCTIONAL DESCRIP'VION 

The functions of the Data Plotting and Fitting Utility (DPFU) are to 
produce Calcomp plots of selected variables from the Measurements Data 
Base and optionally compute fits to the variables. Plots may be 
stacked or overlayed as desired, To make the DPFU flexible and 
convenient to use, most of the important plot parameters are 
determined by options specified in the NAMELIST input file. Among the 
options is the ability to chouse the variables that are to be plotted. 
While there are few inherent restrictions on which variables can be 
plotted against each other, some typical selections include 

0 predicted measurement vs. time 
0 observed measurement vs. time 

o residuals vs. time 

o residuals vs. spacecraft orbit positon 

0 residuals vs. scanner or bolometer temperature 

0 residuals vs. horizon crossing latitude 
0 triggering heights vs. horizon crossings latitude 

Other plot parameters that are set by NAMELIST options are plot 
scales, axis labels, plot title, and number of plots to be stacked or 
overlayed. In addition, the plot can be drawn with or without a data- 
fitting line. 

The following data fitting procedures are available: 

1. Polynomial fit 

2. Finite Fourier Series Fit 
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Also, the system allows a combination of polynomial and finite Fourier 
series terms to be fit. 





While typical use of the fitting procedures is to plot the fitting 
line, there is an option to suppress the plot and just print the 
coefficients of the fit in the Summary Report. It is also possible bo 
plot a "fib" that is input from the NAMELIST and nob computed from the 
selected data span. 

By default, the DPFU computes the mean of each variable and displays 
it on the plot. This feature can be suppressed if desired. 

If the option to stack or overlay a series of plots is selected, each 
plot is defined by a separate NAMELIST file. The DPFU Interates 
through the NAMELIST files until all of the plots are drawn. 





7.2 


INPUTS AND OUTPUTS 


7.2.1 NAMELIST Input 


The Data Plotting and Fitting Utility reads the following parameters 
from the NAMELIST input file, &DPFUIN; 


lOPT 


N PLOTS 
IFIT 


Plot option flag: 
lOPT = 0 
s 1 
= 2 


No Plot 
Standi^rd Plot 
Serial Stacked Plot 


= 3 Parallel Stacked Plot 

s H Overlay Plot 

Number of plots (for lOPT = 2,3,4) 

Data-fitting option: 

IFIT =0 no fitting 

= 1 perform fib and plot fit 

= 2 perform fit 

= 3 use input fit coefficients 


NAVG 
ISIFT 
ICONST \ 
IDEG(5) / 
ISIN(7) ? 
ICOS(7) ) 
PERIOD 
CCONST \ 
CDEG(5) I 
CSIN(7) ( 
CC0S(7) 
XPOINT 


) 


Number of points to be used in data averaging 
Data sift factor for plot and data-fitting procedure 
Arrays indicating selection of basis functions for 
data-fitting procedure. 


Basic period for finite Fourier series fit 

Input fit coefficients for optional plot of "fit" for 
NAMELIST input coefficients. 

Cut-off point on x-axis for series plot (lOPT = 2) 
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IXVAR I X 
lYVAR ) 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 


and Y axis variable indicators, where 

Observed Earth Width 
Observed Earth Phase 
Observed Earth-in angle 
Observed Earth-out angle 
Predicted Earth width 
Predicted Earth Phase 
Predicted Earth-in angle 
Predicted Earth-out angle 
Residual Earth Width error 
Residual Earth Phase error 
Residual Earth-in angle error 
Residual Earth-out angle error 
Residual Pitch error 
Residual Roll error 
Earth-in triggering height 
Earth-out triggering height 
Bolometer temperature 
Scanner temperature 
Signal Status 
Sensor Status 

Predicted Earth-in horizon crossing latitude 

Predicted Earth-out horizon crossing latitude 

Reference Pitch 

Ri^ference Roll 

Reference Yaw 

Spacecraft altitude 

Subsatellite latitude 

Subsatellite longitude 
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29 

Satellite orbit phase angle relative to 
ascending node 


30 

Time 

ISCAN 


Scanner number (1 or 2) 

MAXN 


Maximum number of points to be plotted 

ICHAR 


Plot character, If ICHAR is not specified then a solid 
line is drawn. 

XLEN 


Length of x-axis in inches 

YLEN 


Length of y-axis in inches 

XMIN 


Minimum x-axis value 

XMAX 


Maximum x-axis value 

YMIN 


Minimum y-axis value 

YMAX 


Maximum y-axis value 

TSTEP 


Time step for tic marks on x-axis when x-variable is 
TIME (in seconds) 

TSCALE 


Scale on x-axis when x-variable is TIME (in sec/inch) 

TITLE 


Plot title 

XTITLE 


x-axis title 

YTITLE 


y-axis title 


7.2.2 Dataset Input 

The dataset input to the Data Plotting and Fitting Utility is from the 
Measurements Data Base. For a description of its format and content, 
see Section 2.5. 

7.2.3 Plot Output 

The Calcomp plot utilities write to an output file which is stored on 
tape. The tape is then carried to the Calcomp plot facilities to 
obtain the hardcopy plot output. Figures 7-1 through 7-4 are 
representative of the kinds of plots generated by the DPFU. 
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Figure 7-3 Parallel Stacked Plot 





'1.2. l\ 


Printout 


Summary Report 

The DPFU outputs a summary of the plotting and fitting process. 
Summary Report contains the following information: 

0 NAMELIST parameters 

0 System-genehated plot parameters such as plot scales 
0 Coefficients of the fit 


The 
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SYSTEM STRUCTURE 


The baseline diagram of the DPFU is given in Figure 7-5. The MAIN 
routine is the driver of the system. If multiple plots Ci.e,, 
parallel or serial stacked plots, overlays) are requested, the MAIN 
routine fully processes the first plot, then repeats the entire 
procedure for each of the remaining plots. Since each plot is 
produced from its own NAMELIST file, certain plot characteristics 
(such as axis labels) can be set individually. 

The dataset input file is processed sequentially so that one point is 
plotted for each input record. If data-fittirig is requested, the MAIN 
routine calls DFIT after each record is input to accumulate the sums 
needed to compute the least-squares fit. After all the data has been 
processed, MAIN calls SOLVE to calculate the coefficients for the fit, 
and then calls PLTCRV to plot tne curve. 

The DPFU uses several standard Calcomp utility subroutines that have 
been omitted from the system description for simplicity. 
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SUBROUTINE DESCRIPTIONS 


Each subroutine of the DPFU is described below. 
MAIN 


Purpose; MAIN serves as the driver of the DPFU. 

Procedure; 

(1) Call NAMLST to read in the plot parameters. 

(2) Call INIT to set global parameters and set the origin. 

(3) Call GRAPH to set scales, draw axes, and write labels. 

(i|) Read a data record. 

(5) Call PLTDAT to plot the data. 

(6) If IFIT = 1 or 2, call DFIT to accumulate sums. 

(7) Repeat steps 4-6 fir each data record. 

(8) If IFIT = 1 or 2, call SOLVE to compute the coefficients 
for the least-squares fit. 

(9) If IFIT = 1 or 3, call PLTCRV to plot the data-fitting 
curve. 

(10) Repeat steps 1-9 for each plot. 

Input Parameters; None 


Output Parameters; None 



NAMLST 


Purpose; NAMLST reads the NAMELIST input file &DPFU1N. 
Procedure: 

(1) Perform NAMELIST read. 

(2) Test for errors. 

Input Parameters: None 

Output Parameters; 

lERR Error flag: 

lERR =0 no error 

s 1 read error 

The remaining output parameters are listed by common block 
/CTRL/ 

lOPT Plot option flag 

NPLOTS Number of plots 

IFIT Data fitting option 

ISIFT Step size for data fitting procedure 

XPOINT Cut-off point on x-axis for series plot 

IXVAR x-variable selection 

lYVAR y-variable selection 

ISCAN Scanner number 

MAXN Maximum number of points to be plotted 


/GRAPH/ 


ICHAR Plot character 

XLEN x-axis length 

YLEN y-axis length 

XMIN Minimum x-value 

XMAX Maximum x-value 

YMIN Minimum y-value 
YMAX Maximum y-value 

TSTEP Time increment for x-axis tic marks 

/LABELS/ 

TITLE Graph title 

XTITLE x-axis title 

YTITLE y-axis title 

/BASIS/ 

ICONST I 

IDEG ( Basis function selection flags, 
ISIN ( 

ICOS ) 


INIT 


Purpose; INIT initializes global parameters and determines the 
position of the plot origin. 


Procedure: 

(1) Initialize parameters. 

(2) Set origin: 

YO s YO + YMOVE 

where YMOVE is based on YLEN. (XO is fixed for all plots) 

Input Parameters: 

YO Current plot origin y-coordinate 

YLEN Length of y-axis 

/BASIS/ 

ICONS! \ 

IDEG / Selection flags for basis functions 
ISIN ? 

ICOS 1 

Output Parameters: 

HEIGHT Height of plot in inches 

WIDTH Width of plot in inches 

/DATFIT/ 

IMAP Array mapping basis function selection flags to basis 

functions 

NFIT Order of the fit I 


GRAPH 


Purpose: GRAPH sets up the "frame'' of the current plot: it determines 

the scale, draws the axes, and writes labels. 


Procedure: 

(1) Call SCALE tc determine plot scale. 

(2) Call DRAWAX to draw the axes and tic marks. 

(3) Call TICMRK to label the tic marks. 

(4) Write plot title, axis labels. 

Input Parameters; 


X0,Y0 

Plot origin. 

XMIN 

Minimum x-value 

XMAX 

Maximum x-value 

YIN 

Minimum y-value 

YMAX 

Maximum y-value 

XLEN 

Length of x-axis 

YLEN 

Length of y-axis 

TITLE 

Graph title 

XTITLE 

x-axis label 

YTITLE 

y-axis label 


inches 

inches 


Output Parameters; 

SCALE Scale factor 

XINC X-axis scale increment 

YINC Y— axis scale increment 
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SCALE 


Purpose: SCALE computes the plot scale. 

Procedure: 

(1) Compute Scales: 

XSCALE = WIDTH/(XMAX - XMIN) 

YSCALE = HEIGHT/(YMAX - YMIN) 

(2) Compute number and value of tic marks. 

Input Parameters; 


HEIGHT 

Height of plot 



WIDTH 

Width of plot 



XMIN 

Minimum x-value 



XMAX 

Maximum x-value 



YMIN 

Minimum y-value 



YMAX 

Maximum y-value 



XLEN 

Length of x-axis 

in 

inches 

YLEN 

Length of y-axis 

in 

inches 


Output Parameters; 

XSCALE X-axis scale factor 

YSCALE Y-axis scale factor 

NXTIC Number of tic marks for x-axis 

NYTIC Number of tic marks for y-axis 

XTICS Tic mark values for x-axis 

YTICS Tic mark values for y-axis 

XDIST Distance between tic marks on x-axis 

YDIST Distance between tic marks on y-axis 
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DRAWAX 


Purpose; Draw the x and y axes and tic marks. 
Procedure: 


(1) 

Draw x-axis 


(2) 

Place tic marks on 

x-axis 

(3) 

Draw y-axis 


(4) 

Place tic marks on 

y-axis 


Input Parameters; 


XO, YO 

Plot origin 


XLEN 

Length of x-axis in 

inches 

YLEN 

Length of y-axis In 

inches 

NXTIC 

Number of tic marks 

on x-axis 

NYTIC 

Number of tic marks 

on y-axis 

XDIST 

Distance between tic 

: marks on x-axis 

YDIST 

Distance between tic 

marks on y-axis 


Output Parameters: None 


TICMRK 


Purpose: TICMRK labels the tic marks on the x and y axes. 

Procedure; 

For each axis, perform the following; 

(1) Based on spacing, determine which tic marks should be 
labeled. 

(2) Use XTICS, YTICS to label the selected tic marks. 


Input Parameters; 


XO, YO 

Plot origin 


NXTICS 

Number of x-axis 

tic marks 

NYTICS 

Number of y-axis 

tic markh 

XDIST 

Distance between 

tic marks on 

YDIST 

Distance between 

tic marks on 

XTICS 

Tic mark values 

on x-axis 

YTICS 

Tic mark values 

on y-axis 


x-axis 

y-axis 


Output Parameters; 


None 


PLTDAT 


Purpose; PLTDAT plots a single data point. 

Procedure : 

(1) Call FNDPOS to determine plot coordinates corresponding to 
the data point. 

(2) If ICHARsblank, leave pen down and move to new data point. 

(3) If ICHAR^blank, raise pen, move to new data point, draw plot 
character ICHAR. 

Input Parameters; 


XDAT 

X-value 


YDAT 

Y-value 


XSCALE 

X-axis scale 

factor 

YSCALE 

Y-axis scale 

factor 

HEIGHT 

Plot height 


WIDTH 

Plot width 


ICHAR 

Plot character 


Output Parameters; 


NONE 
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FNDPOS 


Purpose; FNPPOS computes the x and y coordinates of a given data 
point. 

Procedure; 

(1) Compute X-coordinate 

(2) Compute Y-coordlnate 

Input Parameters; 


XMIN 

X-axis minimum value 

XMAX 

X-axis maximum value 

YMIN 

Y-axis minimum value 

YMAX 

Y-axis maximum value 

XSCALE 

X-axis scale factor 

YSCALE 

Y-axis scale factor 

XDAT 

X-value 

YDAT 

Y~value 


Output Parameters; 

X X-coordinate 

Y Y-coordinate 


Purpose; DFIT aocumulates the sums needed for the least-squares fit 
using the selected data points. 


Procedure; 

(1) Accumulate right-hand side vector; 

B(K> = B(K) + F(K,X)»Y FOR K=1, NFIT 

(2) Accumulate least-squares matrix; 

A(J,K> = A(J,K) + F(J.X) « F(K,X) 
for J=1, NFIT; Ksl, NFIT 

Input Parameters; 

NFIT Order of the fit 

X X-value 

Y Y-value 

Output Parameters; 

B Right hand aici« of least squares system 

A Least-square matrix 


F 

Purpose; F is a function that evaluates a selected basis function at 
a given point X. 


Procedure; 

(1) Determine which basis function to evaluate based on input 
parameter IBASIS, 

(2) Evaluate the selected function. 

Input Parameters: 


IBASIS Index indicating which basis function is to be evaluated. 
X Point to be evaluated. 

The remaining input parameters are from common block /DATFIT/; 

IMAP Array mapping IBASIS index to selected basis function. 
NFIT Order of the fit. 


Output Parameters; 

F Function value at point X. 
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SOLVE 


Purpose: SOLVE computes the solution of the linear system of 

equations, AX :: B. 


Procedure; 

(1) Compute the LU decomposition of the matrix A using the 
Cholesky decomposition method. 

(2) Back-solve to compute X. 

Input Parameters: 

A Matrix of coefficients 

B Right hand side 

N Order of the system 

NDIM Actual row dimension of A 

Output Parameters: 

X Solution ’.factor 

lERR Error flag: 

lERR = 0 no error 

s 1 A is singular 
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. .... V.,., 


PLTCRV 


Purpose: PLTCRV plots the dsta-mtlnc curve 


Procedure; 


( 1 ) 

( 2 ) 

(3) 

('O 


Evaluate least-squares fit at selected point. 

Cell FNDPOS to Obtain X and Y coordinates of 'the data point. 
Move pen to computed position. 

Repeat (1)-(3) for each selected point. 


Input Parameters; 


c 

Array of computed least 

XMIN 

Minimum X-value 

XMAX 

Maximum X-value 

YMIN 

Minimum Y-value 

YMAX 

Maximum Y-value 

XSCALE 

X-axis scale factor 

YSCALE 

Y-axis scale factor 

XINC 

Step size for curve 


-squares coefficients 


Output Parameters; 
NONE 
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SECTION 8 - SCAN PATH PLOTTER 


8.1 FUNCTIONAL DESCRIPTION 

The scan path plotter draws the ground path of a scanner field-of- 
view on a perspective view of the Earth grid. This is used to show 
the positions on the Earth viewed by the scanner. The Earth grids 
generated by this utility can be overlaid on GOES infrared photographs 
to show the meteorological conditions viewed by the scanner. 

The following steps are performed in the computation of the scanner 
ground track. 

Step 1 


The perspective view of a spherical latitude-longitude grid is drawn, 
and parameters are initialized for the drawing of general curves on 
this grid. 

Step 2 

The input attitude and spacecraft position information are used to 
compute the scanner spin axis vector and expected Earth width. 

Step 3 

For each scanner rotation angle in the sweep across the Earth, the 
latitude and longitude of intersection of the scanner line-of-sight 
with the Earth is computed and stored in an array. This intersection 
is computed for very amall steps in the scanner rotation angle near 
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the horizons, where the surface intersection position changes rapidly, 
and for larger steps in the scanner rotation angle near the middle of 
the scan. 

Step ^ 


The curve described by the arrays of latitude and longitude 
intersection points is plotted on the Earth grid. 

Step 5 

Steps 3 and 4 are repeated for a second scan cone angle. The two scan 
cone angles are usually selected to bracket the sensor field-of-view. 

Step 6 

Lines are drawn connecting the scan cones at selected positions. 
These positions may be specified by the following: 

0 Degrees of rotation angle from the Earth center direction 
0 Degrees of rotation angle after the Earth-in crossing (specifically 
the average of the Earth-in crossings for the two cones) 

0 Degrees of rotation angle before the Earth-out crossing 
(specifically, the average of the Earth-out crossings for the two 
cones) 

Step 7 

The subsatellite position is marked on the Earth Grid. 
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step 8 


Steps 2 through 7 are repeated for additional scan path plots. 


8.2 


INPUTS AND OUTPUTS 


8.2.1 NAMELIST Inputs 

The scan path plot generation is controlled through two different 
NAMELISTS which are read from the same input file. The first NAMELIST 
which is read, &EARTHG, defines the sphere grid upon which the scan 
paths are to be plotted. The next NAMELIST which is read, &IRSCAN, 
defines the parameters for drawing the ground track of a pair of scan 
cones which bracket the sensor field-of-view. As many of the &IRSCAN 
NAMELISTS as desired may be stacked in the NAMELIST input dataset so 
that several scan paths may be drawn on one Earth grid. A control 
parameter in the &IRSCAN NAMELIST indicates whether a new Earth Grid 
should be drawn. 

Thesf two NAMELISTS are described in the following pages. 
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&EARTHG NAMELIST: 


ALFO 

DELO 

RADIUS 

RVIEW 

RERTHG 

INTVL 

ITVPOL 
L AT POL 


East longitude of the center of the view of the Earth 
grid, from 0 to 360 degrees (West longitudes can be 
converted to East lonitude by subtracting from 360). 
Latitude of the center of the view of the Earth grid, 
from -90 to 90 degrees (South latitudes are given by 
negative values). 

Radius of Earth grid drawing in inches. 

Distance at which the Earth grid is viewed in kilometers 
for perspective effects. 

Radius of Earth grid in kilometers for the computation, 
of perspective effects (Default value is 6378.14). 
Interval of latitude and longitude lines on the Earth 
grid in degrees. 

Interval of longitude lines near the poles in degrees.' 
North and South latitude past which ITVPOL applies. 
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&IRSCAN NAMELIST 


Main Control Parameters 


ISCAN Scan Drawing Option 

= 0, Do not draw scan path 
= 1, Draw pair of scan paths 
ISUB Subsatellite point marking option 

= 0, Do not mark the subsatellite point 
= 1, Draw mark to indicate the subsatellite point 
NEWGRD New Earth grid selection option. 

= 0, Do not draw a new Earth Grid, The next 
&IRSCAN NAMELIST which is read from the input 
dataset will define a new scan path to be 
drawn on the same grid. 

= 1, Draw a new Earth grid. The next &EARTHG 
NAMELIST which is read from the input dataset 
will define a new grid to be plotted. 


Scan Cone Parameters 


SC0NE(2) Pair of scan cone angles for which ground paths 

will be drawn (these are usually chosen to bracket 
the sensor field-of-view) 

IMARKSC200) Indicator for positions on the scan path which are 
to be marked by lines between the two scan cone 
paths . 

= 0, No mark 

= 1, Draw mark at SMARKS(I) degrees of scan 
rotation angle from the middle of the Earth 
width . 
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s 2, Draw mark at SMARKS(I) degrees after the 
Earth-in crossing direction, 
s 3, Draw mark at SMARKS(I) degrees befor? the 
Earth-out crossing. 

SHARKS (200) Scan rotations for marks as defined by the variable 

IMARKS. 

Scan Axis Definition Parameters 

lATT Indicator of the source of the scan axis attitude 

= 0, Input the scan axis right ascension and 
declination, using the NAMELIST parameters 
RAAXI and DECAXI 

= 1 , Specify the scan axis by its azimuth relative 
to the flight path and its elevation relative 
to the local horizontal plane, using the 
NAMELIST parameters AZAXI and ELAXI. 

= 2, Input the spacecraft pitch, roll, yaw attitude 
and the scanner mounting alignment angles, 
using the NAMELIST parameters PITCH, ROLL, 
YAW, AZI, TILT. 

Note: For lATT = x or 2 the orbit parameters EYE and 
RANODE are used to determine the flight 
direction so that these orientations relative 
to the orbit frame are defined in GCI 
coordinates. 

RAAXI Right ascension of the scan axis in degrees 

DECAXI Declination of the scan axis in degrees 

AZAXI Azimuth angle of the scan axis, defined as a right 

hand rotation about the nadir vector from the 
flight direction to the scan axis in degrees 
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ELAXI 

PITCH 

ROLL 

YAV/ 

AZI 

TILT 


Elevatioii of the scan axis above or below the local 
horizontal plane, in degrees 
Spacecraft pitch attitude 
Spacecraft roll attitude 
Spacecraft yaw attitude 
Scanner azimuth alignment 
Scanner tilt alignment 


Spacecraft Position Parameters 


IPOS 


TIME 

RASAT 

DECSAT 

DISTF, 

SATLAT 

SATLON 

SATALT 

EPOCH 


Indicator of the source of the spacecraft position 
relative to the Earth grid. 

= 0, Input the spacecraft right ascension and 
declination distance from the geocenter, and 
G[1T time, using the NAMELIST parameters RASAT, 
DECSAT, DISTE, and TIME. 

z 1 , Input the spacecraft subsatellite latitude and 
longitude and altitude above the surface, 
using the NAMELIST parameters SATLAT, SATLON, 
SATALT. 

= 2, Input the TIME and spacecraft orbital 
elements, using the NAMELIST parameters TIME, 
A, E, EYE, WO, RANODE, EMO, and EPOCH. 

Time in format YYMMDD.HHMMSS 

Spacecraft right ascension position 

Spacecraft declination position 

Distance from the spacecraft to the Earth center 

Subsatellite latitude 

Subsatellite longitude 

Satellite altitude above the surface 

Epoch time for orbital elements in format 

YYMMDD.HHMMSS 


8-8 


A 

e; 

EYE 

WO 

RANODE 

EMO 


Orbit semi-major axis (km) 

Orbit eccentricity 

Orbit inclination 

Orbit argument of perigee 

Orbit right ascension of ascending node 

Orbit mean anomoly 


Earth Parameters 


REARTH 

FLATC 


Radius of the Earth in kilometers 
Earth oblateness coefficient 
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SYSTEM STRUCTURE 


Figure 8-1 shows the baseline diagram for the Scan Path Plotter. 
SCMAIN serves as the driver for the system. It reads the input 

NAMELISTS and controls the program flow. Figure 8-2 shows the data 
flow. 
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Figure 8-2 Scan Path Plotter Flow Diagram 
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8.4 


subrojtine descriptions 


The purpose, procedure, inputs, and outputs of the scan path plotter 

subroutines are described as follows; 

SC MAIN 

Purpose: SCMAIN serves as the driver of the scan path plotter. 

Procedures ; 

(1) Read and write NAMELIST input dataset card images. 

(2) Read the NAMELIST 4EARTHG and call SPHGRD to draw the Earth 
grid. 

(3) Read the NAMELIST &IRSCAN. 

(4) Call TC0N40 and RAGREN to compute the right ascension of 
Greenwich, England, 

(5) Get the satellite position vector in GCI coordinates by 
calling VEC, or ORBGEN according to the spacecraft position 
input option IPOS, 

(6) Get the scan a^is vector in GCI coordinates by calling 
EULZYX, MATMPY, or VEC according to the scan axis attitude 
input option lATT. 

(7) Call SCAN to step the scan across the Earth and compute the 
ground brack of the scanner. 

(8) Call SPHPLT to draw the scan cone on the Earth grid. 

(9) Mark the specified positions on the scan cone. 

(10) If required, plot the subsatellite point marker. 

(11) If NEWGRD = 0, return to step 3 
If NEWGRD = 1, return to step 4 
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stop when an end-of-file is reached on the NAMELIST input 
dataset . 


Input Parameters: 


None 


Output Parameters: None 


Purpose: To compute the Earth in crossing, Earth out crossing and 

ground track of a scan across the Earth. 


Procedure: 

(1) Get an estimate for the Earth width, using the spherical 
Earth equation. 

(2) Choose a scan rotation angle just before the horizon entry 
position. 

(3) Determine the scanner line-of-sight vector for this scan 
rotation angle by calling VPHASE. 

(il) Determine if the line of sight vector intersects the Earth 
and, if it does, determine the position of the nearest 
intersection point, by calling ELIPSX. 

(5) Call RADECM to obtain the latitude and longitude of the 
Earth intersection points and store in the output arrays. To 
keep the array usage down, do not store points within a 
fixed proximity to previous points. 

(6) Increment the scan rotation angle and return to step (3) 
until the line-of-sight leaves the Earth. Th^ ■ /-’crement is 
very small near the horizon crossings and .i,irg<j. near the 
middle of the scan. 

(7) Call RADECM to obtain the latitude and longitude of the last 
Earth intersection point and store in the output _.*rays. 

Input Parameters: 


VSCANR 

GAMMA 

HEARTH 


Scan axis unit vector 
Scan cone radii's 
Earth equatorial radius 


FLATC 

VERSA! 

ILEVEL 


Earth oblateness coefficient 

Earth to satellite position vector in kilometers 
Debug printout level 


Output Parameters: 

PHIIN 

PHIOUT 

SCNLAT(N) 

SCNLON(N) 

N 


Earth in rotation angle 

Earth out rotation angle 

Array of latitudes of the ground track 

Array of longitudes of the ground track 

Number of points in ground track arrays 


ELIPSX 




Purpose: To determine if a particular line of sight intersects an 

Earth ellipsoid, and if it does, to provide the position 
vector of the first intersection point. 

Procedure: The formulas for this evaluation are discussed in 

Reference 4, section 3.6.1. 

(1) Determine if the line of sight intersects the ellipsoid. If 
not, return flag values. 

(2) Determine the distance to the nearest point of intersection. 

(3) If this distance is negative, the line-of-sight points away 
from the Earth or the viewer is inside the Earth and flag 
values are returned. 

(4) Use the distance to the nearest p^int to obtain the position 
vector of the intersection. 

Input Parameters: 

VFOV 
VERSAT 
REARTH 
FLATC 

Output Parameters: 

VHIT(3) Position vector of ellipse Intersection in kilometers 

lERR Error indicator 

= 0, No error 

s 1, line does not hit Earth 

= 2, viewer inside Earth 

= 3, viewer looking away from Earth 

(If lERR .NE. 0, VHIT is set to (2.0, 2.0, a.O)"^) 


Line of sight unit vector 

Vector from the Earth to the satellite 

Earth equatorial radius 

Earth flattening coefficient 


8-17 







REFERENCES 


1 . LAHDSAT-D Conical Scanner Evaluation Plan . Stephen Bilanow, Lily 
Chen, General Software Corporation GSC-TR8203, March 1982, 

2. Horizon Radiance Data Base Format (Revised) . S.P, Singhal, 
Computer Science Corporation 861-1-4260^1-09, Feb, 4, 1982. 

3. LANSAT-D Data Format Control Book . Volume 2, General Electric SVS- 
10123, July 1981, 

4. LANDSAT-D Conical Scanner Mathematical Modeling . Stephen Bilanow, 
Lily Chen, Igor Kulikov, General Software Corporation GSC-TR8202, 
February 1982. 

5. The Horizon Radiance Modeling Utility System Description and 
User's Guide (Premiminary Draft). Dr, W, Nutt, Computer Sciences 
Corporation CSC/SD-78/6032, March 6, 1978. 


R-1 


t ^ 


